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\ (54) TUtei METHOD AND SYSTEM FOR DISEASE DBTECnON USING MARKER COMBINATIONS 




(57) Al>stract: Hie present invention relates to methods and system for the di^osis diseases or conditions. In a particniar aspect, 
a disclosed method for determining a panel includes calculating a panel response for each patient in a set of diseased patients and in 
a set of non-diseased patients. The panel response is a function of die value of each of a plurality of markers In a panel of marfcere. 
The method also includes calculating a value for an objective function. The objective function is indicative of the effectiveness of 
the panel. The steps of calculating a panel response for each patient and calculating a value for an objecliVB fimction are iterated by 
varying at least one of the parameters relating to the panel response function and a sense of each marioer to facilitate optimization 
of the objective function. The objective function may be a measure of an overlap of panel responses of diseased patients and panel 
responses of non-diseased patients. The contribution of each marker to the objective function may be detennined, and the panel size 
may be reduced by removing the poorest markers. Thus, an optimum panel of markers and an optimal panel response function for 
the diagnosis of a disease or condition may be determined. 
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METHOD AND SYSTEM FOR DISEASE DETECTION 
USING MARKER COMBINATIONS 

FIELD OF THE INVENTION 

[OOOIJ The present inveatioii relates to the identification and use of diagnostic markers 
for various diseases or conditions. More particularly, the invention relates to methods and 
systems for identifying and utilizing panel of markeis far detectloii of one or more 
particular diseases or conditions. 

BACKGROUND OF THE INVENTION 

[0002] The background of the invoition is provided to aid the reader in understanding 
the invention and is not admitted to describe or constiftite prior art to the present inventioda. 

[0003] The clinical presentation of certain diseases can often be strikingly similar, 
even though the underlying diseases, and the qrpropriate treatments to be given to one 
suffering from the various diseases, can be com5)letely distinct. For example, subjects may 
present in an urgent care facility exhibiting a deceptively simple constellation of apparent 
synq>toms {e.g., fever, shortness of breath, dizziness, headache) that maybe characteristic 
of a variety of unrelated conditions. Diagnostic methods often involve the comparison of 
symptoms and/or diagnostic test results known to be associated witii one or more diseases 
that exhibit a similar clinical presentation to the synq)toms and/or diagnostic results 
exhibited by the subject, in order to identify the underlying disease or condition present in 
the subject. 

[OOOi] The acuteness or severity of the sym|>toms oftai dictates how rapidly a 
diagnosis must be estabUshed and treatment initiated. For exan^le, immediate diagnosis 
and care of a patient experiencing a variety of acute conditions can be critical. See, e.g., 
HaniSj iittsfc Fam. Physician 31: 802-06 (2002) (asthma); Goldhaber, Eur. Respir. J. 
Sm>l. 35: 22s-27s (2002) (puhnonaiy emboUsm); Lundergan et ai. Am. Heart J. 144: 
456-62 (2002) (myocardial inferction). However, even in cases whore the apparent 
symptoms appear relatively stable, rapid diagnosis, and tiie r^d initiation of treatment, 
can provide both relief from immediate discomfort and advantageous in^rovement in 
prognosis. 
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10005] Recently, woricers seeking to provide rapid diagnostic methods for various 
diseases or conditions have sought to identify "markers" for diseases; that is, molecules 
that are present in a sample obtained from a subject suffering fiom a disease of interest in 
an amount that dijHers fiom the amount preseait in a sample fixwn a "normal," non-diseased 
object. 

[0006] Diagnoses of many diseases or conditions, such as cardiovascular disease and 
stroke^ fiir example, arepeidaraied by measurement of the levds of particular maikos in a 
patieiit. Often, however, a single marker is generally incapable of providing clinical utility 
because its value does not provide a means of confidently distinguishing between a 
diseased patient and a non-diseased patiait. 

[0007] As an example. Figure 1 illustrates liiat the levels of a particular marker 
eapressed in a diseased and a non-diseased population. As shown in the figure, the masker 
levels in these two populations may be distributed over broad ranges in a distribution 
pattern. Although the diseased population in fbis example generally may exhibit higher or 
lower levels for the marker than the non-diseased population, substantial portions of each 
population M within a region of overlapping values. Thus, definitive or confident 
diagnosis of a disease or a condition based on the measuremeait of tins single marker may 
be impossible. Traditionally one chooses a cutoff value in the overlap regioiL The cutoff 
is chosen to optimize the munber of &lse positive versus the number of fake negatives. In 
practice physicians often treat a patient based on where fhey fall relative to the cutoff. 
They often do not consider how close 1he patient is to the cutoff. 

[0008] The effectiveness of a test having such an overlap is often expressed using a 
ROC (Receiver Operating Characteristic) curve. Other measures, such as positive 
predictive value (PPV) and negative predictive value (NPV) may also be used as a 
measure of the effectiveness of the test ROC curves are well known to those skilled in 
Ibe art. Thus, the details pertaining to ROC curves are beyond the scope of this docume».t, 
however there is a brief descx^ition below. Further, refereoice may be made to Zweig, 
MH. & Canqjbell, C.C., Clin Chem 32, 561-577 (1993) and Hcndrson, A.R., Ann. Clin. 
Biochem3£l, 521-539 (1993). 

[0009] Figure 3 illustrates an example of a ROC curve for the marker level 
distributions of Figure 1. The ROC curve shows the trade off betweai the sensitivity and 
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Specificity of a marker. The set^itivity is a measure of the ability of the marker to detect 
the disease, and the specificity is a measure of the ability of the marker to detect the 
absence of the disease. The horizontal ads of the ROC curve represents (1- specificity), 
which increases with the rate of false positives. The vertical axis of the curve represents 
sensitivity, which increases with the laite of true positives. Thus, for a particular cutoff 
selected, the vahies of specificity and sensitivity may be determined. The right hand end of 
the curve is the Tninimnnri cutoff^ the lefi hand end of the curve is the maofi'mnm cutofif. As 
file cutoff is changed to increase specificity, sensitivity usually is reduced and vice versa. 
The area under the ROC curve is a measure of the utility of tiie measured maikcr level in 
the correct identification of one or more diseases or conditions. Thus, the area under the 
ROC curve can be used to detennme the eSectiveness of the test. Note the area is 
independent of the cutoff value. 

[OOIG] Panels of multiple maakers may improve the likelihood of an accurate 
diagnosis. The multq)l6 marker "panel" for a particular disease is preferably selected such 
that a particular "profile" of marker levels is specific fijr that disease and capable of clearly 
distinguishing disease from non-disease. However, methods fox identifying such panels, 
and tiie particular "profiles" that provide clinical utility, are typically empirical m nature, 
relying on trial-and-enxtr. Furthermore, becaure tiie computational complexity involved in 
identifymg suitable diagnostic tiiresholds and/or profiles mcreases as the number of 
markers in a potential panel inoiease, marker panels typically involve only a few markers. 
Searching for an eflBactivc panel from among a large number of markers can become the 
coinputational equivalent of finding a needle in a haystack. For example, often one might 
look for elevation of 4 of 6 markers, or more generally n of m markers, to define a positive 
state. In this example the cutoff values for each marker are chosen, tiien the data analyzed 
to see how effective the test is. This is repeated for dififerent number of elevated markers, 
cutoffs and markers. In this exanq>le, all madcers are treated with equal importance, there 
is no method to adjust the relative importance. 

BRIEF SUMMARY OF THE INVENTION 

[0011] The metiiod disclosed m tiiis document provides a means to systematically find 
the optimal majkeis and panels of markers to distinguish (compare) non-disease fix>m 
disease, and it also optimizes the way in which die marker values are used. A first step to 
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sdmplify the problem of defboing a marker or a panel of markers is defining an 'objective 
function'. An objective function is a scalar fimction, and will represent the effectiveaiess of 
the test for diagnosis of non-disease from disease. For ejtample, raflier than requiring n 
elevated markers to define a positive state and then quantifying the efifectivene^ of this 
algorithm, one can generate a ROC curve fiom the number of elevated maiketrs, and use 
flie area undo: flie ROC curve Cthe ROC curve area") to define the effectiveness of tiie 
test. By using tiie ROC curve area m the effectiveness of the test, the optimization 
problem has been sinqilified. This is because the search Efpace has been reduced since there 
is no need to calculate the efEectiveness associated with each of the m values for n 
elevated markers, hi this example, the number of elevated markere can be thought of as a 
concentration for the ROC curve, but as described above, the selection of the cutoff 
concentration is not required to determine If a test will be effective. Another step to 
simpHjfy the problem of defining a marker or a panel of markers may be to define a 
systematic way to find the best way to use the markers. Without this it is very difficult to 
find the best markers because one needs to distinguish the markers and how to use them. 
A systematic method to find the best way to use the markers is to combine all the values 
into one result, the "panel response". Functional forms of the panel response can be 
selected in a variety of manners described hcmn. Once this is done search routines can be 
employed to find the panel response fimction to maximize or minimize the objective 
function for a set of madcers. 

[0012] The m^od may also inchide a technique for determining the relative 
importance of the markers in the set, and subsequently determine the optimum maikears to 
use, for example, in a panel of n markers. 

[0013] In addition to measured marker levels, other information including a patient' s 
history, sex, age, nice, and oilier factors may also require consideration. In tins regard, 
embodiments of tiie disclosed me^d may accommodate such factors as madcers. 

[0014] SpecificaUy, certain disclosed embodiments ofthfipresesatinveation relate to 
tiie identification and use of diagnostic matkexs for cardiac diseases and stroke and 
cerebral mjury. Generally, the methods and systems described herdn can meet the need in 
the art for the development of an elective panel of markers for the accurate diagnosis of a 
selected disease or condition. More generally, the disclosed metiiods and systems maybe 
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used to develop criteria for distinguishing memibffls of two or more gToig)s fijr whom the 
distribution of certain characteristics are known. 

[00151 ^ ^ fti^t aspect, tiie invention discloses a method of identifying a panel of 
markers for detemmiing a diagnosis of a disease or a condition, or of a progaosis (that is a 
risk of some future outcome). The method inchides calculating a panel response for each 
patient in two sets of patiaits, referred to for conveniaice as "diseased" patients and "non- 
diseased" patients. The panel response is a function of value of each of a plurality of 
matk^ in a panel of maikero. 

[0016] Thetenn"i}anera8U8edha:einreferstoasetofmaik)efs. The panel may 
include any piractical number of maik^ appropriate fin: use with the diagnosis of the 
particular one or more diseases or conditions. 

[0017] Hie toDtn "marker" as used herein refers to proteins, polypeptides, nucleic 
acids, bacteria, viruses, prions, small molecules and the liks^ to be used as targets for 
screening test samples obtained jfrom subjects. "Protems, polypeptides, car small 
molecules" used as markers in the present invention are contemplated to include any 
fragments tiiereof, in particular, immunologically detectable ftagments. "Market'*, as used 
herein, may include d^ved markers as defined below, and may also include such 
characteristics as patient's history, ag6> sex and race, for example. Certain markers are 
also known in the field as "aoalytes". A marker is said to be a scpedfic marker of the 
disease if only the presence or absence of the target disease condition influences its value. 
A marker is said to be a nonspecific marker of the disease if many disease conditions 
influence its value. An example of a specific marker is cardiac Tnl, which, when elevated 
above about I ng/ml is specific to myocardial infarction. An example of a non specific 
marker is C-reactive proteni ("CRP")» which is elevated in conditions that promote the 
inflammatory response. 

[0018] The term "related maricer" as used herein refisrs to one or more ftagments of a 
particular maiker or its biosynthetic parent that may be detected as a surrogate for the 
marker itself or as independent markers. For example, human BNP is derived by 
proteolysis of a 108 amino acid precursor molecule, referred to hereinafter as BNPmos. 
^/bltllre BNP, or "the BNP natriur^c peptide," or "BNP-32'' is a 32 amino acid molecule 
resenting amino acids 77-108 of this precursor, which may be refiarred to as BlSIP77.,o8- 
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The remaining residues 1-76 are refered to hereinafter as BNP1.76. AdditionaUy, related 
markers may be the result of covalent modification of tiie parent marker, for example by 
oxidation of methionine residues, ubiqmtination, cysteinylation, nitrosylation, 

glycosylation, etc, 

[0019] Preferably, the mediods described heremafter utilizB one or more maricers that 
are derived from the subject. The term "subgect-derived maiker" as used herein refers to 
protein, polypeptide, phospholipid, nucleic acid, prion, or small molecule markers that are 
expressed or produced by one or more cells of the subject. The presence, absence, amount, 
or change in amount of one or more maikeis may indicate that a particular disease is 
present, or may indicate that a particular disease is absent Additional markers may be 
used that are derived not from the subject, such as molecules expressed by pathogenic or 
infections organisms that are correlated with a particular disease, race, time since onset, 
sex, Such marirers are preferably protein, polypeptide, phosphol^d, nucleic add, 
prion, or small molecule markers that identify tibe infectious diseases described above. 

[00201 The temi "diagnosis" as used herem refers to methods by which tlie skilled 
artisan can estimate aiid/oT deterrnine whether or not a patient is sufiS^ 
disease or condition. The skilled artisan often makes a diagnosis on the basis of one or 
more diagnostic maricers, the presence, absencei, or amount of which may be indicative of 
the presencei, severity, or absence of the comditioiL hi addition to markers, other tests, 
such as ECG, Echo, and MRI, and other factors, such as patlexit's history, sex, age, and 
race, may also be used in making the diagjK^is. As used heraiD, the term "markers" also 
includes these other tests and other &cftors. While the instant specification describes the 
invention in terms of diagnosis of diseasei, the methods described herein are equally 
appUcable to identifying maricers fin- use m prognosis. A "prognosis" may be determined 
by measuring one or more maricers, the presence or amount of which in a subject (or a 
sample obtained &om Ihe subject) signal a probabiUty that a given course or outcome will 
occur. 

[0021] The tem "panel response" as used herein refers to a scalar fimction or its 
value, which is a fimction of the marker values of the panel. Most generally, the panel 
response is a function of the markra: values (Mi_„) , written asPJ! = ) . In a 
prefeaxed eanbodunent the panel response is a summation over mdicator values (J) of each 
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maiker. The indicator value is generally a fimction of the marker value. This can be 
represented as PR = '^I,{Mi)- JF, , where Jj is a fimctica of the marker valueJW", , If, is a 

weighting coeflBicieflt that scales the indicator function. For definitive purposes, in this 
document it will be assumed that the panel response is scaled such lhat all values are 
betv/eeu 0 and 1, but other mcrements can apply. 

[0022] The tems "diseased" and "non-diseased" as used herein refer to two 
populations that differ in terms of a disease characteristic of interest. For example^ a 
"diseased" population may be a population suffering fiom a stroke, while a. "non-diseased" 
population may be a population not suffering fixim a stroke. Hiese are simply labels, and 
the skilled artisan will understand that one could just as easily label the population 
suffering from a stroke "non-diseased" and the otha: population "diseased" without 
varying from the teachings set foth herein. In flie case of a prognosis, the two populations 
may each be suffering from the same disease, but sei>arBted mto populations on the basis 
of one group having a particular outcome (e.g., death) and the second groiqj not having 
tiiat outcome. The set of diseased patients and set of non-diseased patients may include 
patients whose state, whether diseased or non-diseased, has been confirmed and for whom 
maiker levels are available for one or more markers. 

[0023] The term "marker value" as used herein refers to a numeric value, such as a 
value representing the result of an assay of the marker. For example, the marker value may 
be expressed in xuiits of concentration or number. When the marker represents 
characteristics such as a patient's history, thai the value may be a numeric represoatation, 
or mq)piDg, of the history infiinnatioiL 

[0024] The term "derived maikei'' as used h^iein refais to a value that is a fttnction of 
one or more measured mariceis. For example, derived markers may be related to the 
diiange over a time interval in one or more measured maricer values, may be related to a 
ratio of measured marker values, may be a marker value at a di£Eerent measurement time, 
ctr may be a complex function such as a panel response function. 

[00251 The metiiod fraHier comprises calculating a value for an objective function, the 
objective function being indicative of an efifoctiveness of Ihe panel 
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[0026] The tarn "objective ftmction" as used herein refes to a scalar function or its 
value, which may be a function of tiie plurality of pand resptmses and known disease 
states or diagnoses of a collection of patient isaii5)les. The objective function is a measure 
of tile clinical e£Eectiveness of the test, or the ability to distinguish disease fiom non- 
disease. An example of an objective fimction is the area under the ROC curve. The 
objective fimction may be related to tiie amount of overlsq) between the diseased and non- 
diseased panel response values. The objective ftmction is a scalar value, which is 
indicative of the eflfectivaciess of the panel. The objective fimction may be defined by a 
user as a fimction of various outputs, such as ROC curve features defined below, of the 
panel responses for the groups of patieasts. 

[0027] The method of the first aspect of tiie invention also comptises iterating Hie 
calculating a panel response for each patient and calculating a value for an otrjective 
function by varying at least one of parameters relating to the panel response fimction and a 
seaset of each marker Id facilitate optimization of tiie objective function. 

[0028] "Iterating" may include repeating tiie steps witii variations in the inputs, where 
the variations may be dependant on tiie outputs of tiie previous iteration. "Varying" may 
include tweaking a parameter by either a predetexmined amoimt, an amount dependant on 
an output of the previous iteration or a random amount 

[0029] The term "sense" as used herein refers to the direction of the response of a 
marker with disease state. If a marker value is elevated in diseased patients relative to non- 
diseased pati^ts, then the maiker is said to have a positive sense. If the marker value is 
lower in diseased patients relative to non-diseased patients tiieai tiie marker is said to hove 
a negative sense. If the probability of a finding the marker value near some specific value 
is elevated in diseased patients relative to non-diseased patients, the sense is said to be 
positive, tf the probability of a finding a maiker value near some specific value is reduced 
in diseased patients relative to non-diseased patients, tiie sense is said to be negative. One 
skilled in the art will recognize that it is trivial to invert functions or map the maiker value 
such that a negative s^e marker can be analyzed in the same way as a positive sense 
madKr. Moreover, as discussed above, the labeling of a population as "diseased" or "non- 
diseased" is merely a question of labeling of two populations for the presence or absence 
of a characteristic of interest Throughout this document the maiker sense is described as 
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positive. Has is for conciseness only, all concepts and claims can apply to both negative 
and positive sense markers, and both positive and negative senses are msplkifiy included. 

[0030] The tram '^arametaB" as VKed hscem refcacs to coefficients, powers, etc. of a 
function that may be varied to modify the functional value. For exaa5)le, if the function is 
a ramp fimction, tiie low threshold and the high threshold, may be two parameters that are 
varied. If the ftmction is a Gaussian the width and location may be two parameters that 
are varied. The optimization process will modify one or more of title par^eters of the 
panel response function, which in one ambodnnejit may Include all of the parameters of 
the used indicator fimctions and weighting coefficients. 

[0031] Accordmg to another aspect of the invsition, a system for identifying a panel 
of markers for dia^osis/prognosis of a disease or a condition includes means for 
calculating a panel response for each patient in a set of diseased patients and in a set of 
non-diseased patients. In one embodiment the panel response is a function of a value of 
each of a plurality of markers in a panel of markers. The means for calculating may be a 
cmtral processing unit (CPU), as may be available on a desktop computer, a l^op 
computer, a workstation or a mainframe, for sample. 

[0032] The system further includes means for calculating a vahie for an objective 
function. The objective function is indicative of the effectiveness of the panel. In certain 
embodiments, an objective function may be a measure of overly of panel responses of 
diseased patients and panel res|>Qiise8 of non-diseased patients. 

[003^ Further, the system include nxeans for iteratively activating the ineaiisfi^ 
calculatiiig a panel response and Ihe means for calculating a vahie for an objective 
function, by varying at least one of the folbwing parametears to faciHtate optimization of 
said objective fimction: parameters relating to tiie panel response function and a sense of 
each marker. 

[00341 in another aspect oftiiemvention, a program product includes tnacbine 
readable program code for causmg a machine to perforan certain melhod steps. The 
meihod steps mclude calculating a panel respotose for each patient m a set of diseased 
patients and in a set of non-diseased patients. The panel response is a ftmction of vahie of 
each of a plurality of markers in a panel of markers. 
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[0035] The method steps further include calculating a value fiwr an objective function. 
The objective fimction is indicative of tiie effeotiveoess of tiw panel. Further, the method 
steps include iterating the steps of calculating a panel response for each patient and 
calculating a value for an objective function by varying at least one of iJje following 
parameters to fecilitate optimizatian of said objective function: paramdCTS relating to the 
panel response ftmction and a sense of each marker. 

[0036] In a preferred cxobodimeaat, the program product includes machine readable 
code embedded in a portable metca-. The term *^>ortable meter," as used horein, may 
include any number of devices having the ability to execute coded instructionB. hi a 
further preferred embodiment, the portable meter is a flaorometer. In an altemate 

embodiment, the portable meter is a reflectometer. 

[0037] In a preferred embodiment, the program product inchides machine readable 
code embedded in a computer. In a further prefened embodiment, the computer is a 
portable computer. In another preferred embodiment, iJie computer is adapted to be 
accessed throu^ a network, such as a public netwodc like the Xntemet. 

[0038] In another preferred embodiment, the computer is adapted to be coupled to an 
analyzer. In a further preferred embodiment, the analyzer is an inamimbassay analyzer. In 
an alternate embodiment, the analyzer is a single nucleotide polymorphism detector, hi 
another embodinient, the analyzer is adapted to sort and count sinoilar and diffe^^ 
particles and cells. 

[0039] In a preferred embodiment, the panel response is a function of the value of an 
indicator for each of a plurality of markers in apanel of markers and a weighting 
coefficient for each marker. The indicator is a mapping, for each of the plurality of 
markers, of marker levels. The mapping is according to an indicator function. The 
iterating includes varying at least one of die weighting coef&cimts, parameters relating to 
the indicator function, and a sense of each maiker to &cilitate optimization of the 
objective fimction. 

[0040] The terra, "indicator fimction" as used herein refers to a scalar function or its 
value, which is a function of a maiker value. The mapping is in accordance with a 
indicator fimction. The indicator function may be any fimction providing a value 
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dependent on the marker level. The indicator function may be a aiding of marker values 
into values that may be more closely related to the probabiHty of diseased state at that 
maiker value. The indicator function may be scaled such that aU values are between. 0 and 
1 . In this document it -will be assumed liiat the indicator fimction is sc^ed such that all 
vahifis are between 0 and 1. This scaling does not infhunce the result of the method, 
howevo: in practice it does simplify some fimnulations. For example, to change a positive 
indicator fimction O^IF) to work with a negative sense maiker the negative indicator 
fimctirai (NIF) may be defibaed ^ NIF = 1 - PIF. 

10041] The temi "mappm^ as used herem refers to a relation betwe«i a value in one 
domain to a value in another domain. The m^ing relation may be a on&4o-one 
relationship or a one-to-many relationship. 

[0042] The term "elevation indicator function" as used herein refers to a scalar 
function ttiat has a high and monotonic rate of change between low and high threshold 
values, and a smaller rate of change elsewhere. Examples of this type of function include 
step, ramp, 'S' or sigmoid functions. One skilled in the ait will recognize that there are 
many such fimctions. 

[0043] The term "localization indicator fimction" as used herein refera to a scalar 
fimction that is peaked near some expected value, and decreases when the marker value is 
fiirttier away from the e?q)ected value. Examples of this type of fimction include triangle, 
square, trapezoid, or Gaussian functions. One skilled in the art wiU recognize that ttioe ace 
many such functions. 

[0044] The term "contributioai" as used herem refers to the relative amount that a 
maiker contributes to the objective fimction. The contribution may be related to the 
in^)ortance of a maiker. 

[0045] The tetrn "test" as used herein refers to a method perfimned i;^ch yields an 
output related to a clhsical outcome. A test may comprise values of 1 or moie markets. A 
test may also be a procedure used in the determination of a panel response. Commonly, a 

test is also an immunoassay. 

[0046] In the method the marker values may be combined into one value, the panel 
re^onse. As described above, in a pre^ted embodiment the panel response is lepiesented 
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as PR= 5]/, (Mj ^^^.Choosiiig different fimctiond 

the way a marker is used. For example, when several nonspecific maifcors are used, flieii 
combined elevations of the markers may indicate a diseased state. The iqppiopriate 
indicator functions could be devation indicator functions as defined above. In this 
example, when the maiker value is below a low threshold Ihem there is little or no change 
in the uidicator Amotion with maiker value, and v^3m above a high threshold than again 
there is little or no change in indicator value. Between these thresholds Ihe indicator vahie 
increases or decreases with maiker vahie. One skilled in the art will recognize there are 
many functions that have Ihis propaly. Physically one can associate the Ihresholds with 
the lower and upper end of the ovorlap region ss illustrated in Figure 4. A panel response 
may be a numerical value fixr each patient The range of values of the panel response may 
be set to any desired range. For example^ the values of the panel response may be scaled to 
fall between zero and one. 

[0047] In another embodiment the indicator function is chosen to localize the marker 
value. For example if a certam pattern of maikBr levels is associated with a disease state 
then the indicator function could be a locaUzation indicator fimction as defined above. 
These functions give a hig^ response when the macko: is near the optinial value. One 
skilled in the ait will recognize there are many functions that have Ihis localization 
property. In an exanaple using 1hese fimctions, certain disease states such as unstable 
angina, may be an intermediate disease. A marker such as Tnl is elevated by minor 
necrosis due to ischemia associated with unstable angina, but is elevated still further by 
major necrosis associated with myocardial infarction. Other maricers may be elevated with 
unstable angina, but not elevated with myocardial in&iction. The mdicator function of 
each analyte can be different so panels can consist of markers of bolfa types, as needed in 
the exaax^le above. 

I004S] hi a preferred embodiment the method includes utilizing a search engme to find 
optimal parameters fbr the panel response function. The search engme is able to efficiently 
vary parameters of the panel response until it finds a set that results in a local 
maximization of tiie objective function. Becaxise the objojtive function is a measure of the 
effectiveness of flie test, Ihe optimized panel response may provide an improved 
diagnostic value. 
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10049] In a preferred embodiment the method includes calculating a contribution for 
each marker. In another preferred embodiment the contributions of all markers are ranked, 
and markers with low contribution values may be removed &om the panel. The entire 
process can be repeated with the reduced numbar of markers until the desired panel size 
and perfonnance are achieved, 

[0050] Another embodiment of the invention measures multiple markers from a 
patient and combines the values into a single panel response. The panel response function 
could be determined by the method described above. The panel value would be compared 
to a cutoff value, providing an effective tool to aid in the diagnosis of disease states. 

[0051] In a prefixed embodiment, an objective flmction is a measure of overlap of 
panel responses of diseased patients and panel responses of non-diseased patients. 

[0052] According to a preferred embodiment, the calcukting of a value for an 
objective fimctioii includes generatiiig a receiver operating cbatactedstic (ROC) curve for 
the panel response. The ROC curve is indicative of a sensitivity of the panel response as a 
fenction of one minus a specificity of the panel response. ROC curves are well-known to 
those sikiUed in the art sad are further described below. 

[0053] In various aspects, multiple deteinnination of the marker panels described 
herein can be made^ and a tenqKual change in fhemaricers can be used to rule in or out one 
or more diagnoses or prognoses. For eacample, one or mote markers may be detcnnined at 
an initial time, and again at a second time. In such embodiments, an increase in the 
madcer fiom the initial time to the second time may be diagnostic of a particular disease, 
or indicate a particular prognosis. Likewise, a decrease in the marker from the initial time 
to the second tune may be indicative of a particular disease, or of a particular prognosis. 

[0054] In yet other embodiments, multiple detemmationsofmaiker panels can be 
made, and a temporal change in the marker can be used to monitor the efficacy appropriate 
tileries. In such an embodiment, one migfal tsxpect to see a decrease or an increase in the 
marker(s) over time during the course of effective therapy. 

[00S5] In yet a further aspect, the invention relates to devices for analyzing the marker 
panels described heajein. Such devices preferably contain apturality of discrete^ 
uulepend^tly addressable locations, or "diagnostic zones," each of 'niiich is related to a 
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particular maiker of inteaBst Following reaction of a sample with the devices, a signal is 
gaierated from flie diagnostic zQne(s), which may then be correlated to fbs presence or 
amount of the markets of intaest In preferred embodiments, one or more of the 
diagnostic ames comprise an antibody that binds for detection the particular marker to be 
detected at that particular addressable location. 

[0056] The term "discrete" as used herein refers to ateas of a surface that are non- 
contiguous. That is, two areas are discrete from one another if a border that is not part of 
either area completely surrounds each of the two areas. 

[0057] The tenn "independently addressable" as used herem refers to discrete areas of 
a sur&ce from whidb a specific signal may be obtained. 

[0058] The term "antibody" as used herein refers to a peptide or polypeptide derived 
from, modeled after or substantially encoded by an immimoglobulin gene or 
inmiunoglobulin genes, or fragments tiiereof, enable of specifically binding an antigen or 
epitope, See. e.g. Fundamental Immunology, S"' Edition, W.E. Paul, ed.. Raven Press, 
N.Y. (1993); Wilson (1994; J. ImmunoL Methods 175:267-273; Yannush (1992) J. 
Biochem. Biophys. Methods 25:85-97. The term antibody includes antigen-binding 
portions, i.e., "antigen binding sites," (e,g., fragmmts, subsequences, complementarity 
determining regions (CDRs)) that retain capacity to bind antigen, including (i) a Fab 
fragment, a monovalent fragmrait consisting of the VL, VH, CL and CHI domains; (ii) a 
F(ab')2 fragment, a bivalent fragment comprising two Fab fragments linked by a disulfide 
bridge at the hinge region; (iii) a Fd fragment consisting of the VH and CBl domains; (iv) 
a Fv fragment consisting of the VL and VH domains of a single arm of an antibody, (v) a 
dAb fragment (Ward et al.. (1989) Nature 341:544-546), which consists of a VH domain; 
and (vi) an isolated complementarity determining region (CDR). Single chain antibodies 
are also included by reference in the term "antibody." 

BRIEF DESCRIPTION OF THE DRAWINGS 

[00S9] In the following, the invention will be explained in fiulher detail with reference 
to the drawings, in 'n^ch: 
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[(K)60J Figure 1 is a chart illustrating aa exemplary distribution of levels of a particular 
marker among a set of diseased patients and a set of non-diseased patients; 

[0061] Figure 2 is a chart illustrating an exeo^laiy scatter distribution of levels of a 
particular marker among a set of diseased patients and a set of non-diseased patients; 

[0062] Figure 3 is an exemplary receiver operating characteristic (R(XO curve i^^ 
marker level distributions illustrated in Figure 2; 

[0063] Figure 4 is illustrates Hie chart of Figure 1 with the raarkear values being 

m^ed to an indicator value; 

[0064] Figure S is a chart illtistrating an exen^lary scatter distribution of panel 
responses for the set of diseased patietots and the set of non-diseased patients; 

[0O65] Figure 6 illustrates a ROC curve for the panel response distributions of Figure 
5 with the knee of the ROC curve labeled; 

[0066] Figure 7 illustrates the progression of ROC curves through an optimization 
process; 

[0067] Figure 8 is a chart illustrating the relative contributions of each marker in a 
panel; 

[0068] Figure 9 shows the individual ROC curves and areas for each of 5 markers 
comprising the panel for figures 6 and 16; 

[00691 Figure 10 shows the initial and final ROC curves for an optimization of 38 

markers; 

[0070] Figure 1 1 shows the ranking and relative average contributions of 38 markers 
after SO optimizations; 

[0071] Figure 12 shows the initial and final ROC curves for an optimization of 19 

markers; 

[00721 Figure 13 shows tiie ranking and relative average contributions of 19 markers 
afta: 50 optimizations; 



IS 



wo 2004/058055 



PCT/US2003/041426 



[0073] Figure 14 sSiows the initial and final ROC cairves for an optinuzation of 10 

markers; 

[0074] Figure 15 shows the ranking and relative average contribiriions of 10 markers 
afier 50 qptimizatiQns; 

[0075] Figure 16 shows the initial and final ROC curves for an optimization of 5 

markers; 

[0076] Figure 1 7 shows the ranking and relative average contributions of 5 markers 

after 50 optimizations; 

[0077] Figure 1 8 shows the optimized ROC curves of 6, 3, and 2 measured and 
derived markers and 3 measured markers used to diagnose AMJ; and 

[0078] Figure 1 9 shows title relative contributions of all 6 of the measured and derived 
markers for AME, 

DETAILED DESCRIPTION OF THE INVENTION 

[0079] In accordance witia the present irivention, there are provided methods and 
systems &r the identification and use of a panel of markers for the diagnosis and/or 
prognosis of one or more conditions or diseases, such as cardiovascular diseases and 
strokes, in a subject 

[0080] Method for Defining Panels of Markers 

[0081] la practice, data may be obtained from a group of subjects. The subjects may 
be patients who have been tested for the presence or level of certain markers. Such 
markers are well known to those skilled in the art. A particular set of markers may be 
relevant to a particular condition or disease. The method is not dependent on the actual 
markers. The markers discussed in this document are included only for illustration and are 
not intended to limit the scope of tiie invaitioa Examples of such markers and panels of 
maiieKi are described in pending U.S. patent application Serial Number 10/139,086, 
entitied "DIAGNOSTIC MARKERS OF ACUTE CORONARY SYNDROMES AND 
METHODS OF USE THEREOF," and U.S. Patent Application Serial Number 
10/225,082, entitled "DIAGNOSTIC MARKERS OF STROKE AND CEREBRAL 
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INJURY AND METHODS OF USE THEREOF," each of winch is assigned to the 
assignee of&o present apphcati<m and is incoiporatied herein fay le&T^ In accordance 
wi3i the disclosed embodimcaits of Ihe pcesetA inveaitioii, '*tnadcers" may also include 
derived maikers as defined above, and :^ctois such as a patient's faistoay, sex, age and 
lace, for example. 

[0082] The group of sobjects is divided into at least two sets. The first set includes 
subjects who hove been confirmed as having a disease or, mote genially, being in a first 
condition state. For eExample, this first set of paiiesits maybe fbose that have recently had 
a stroke. The confinnaticaof&is condition state inay be made through more rigorous 
and/or expensive testing. For purposes of this document, it will be assumed that this « 
testing is able to confirm the condition state. Hereinafier, aibjects in this first set will be 
referred to as "diseased". 

[0083] The second set of subj ects are selected fiom those who do not faU within the 
first set. Has set may inchide all remainmg subjects, or only those subjects being in a 
second condition state. Subjects in this second set will hereinaft^ be re£arred to as "tion- 
diseased". Preferably, the first set and the second set each have an approximately equal 
number of subjects. The first and second sets of data are said to be a group of data. 
Multiple groups of data maybe defined by rqi^ding the steps above for diSerent disease 
states, condition states, or any other selection criteria. 

[0084] The data obtained fi-om subjects in these sets includes levels of a plurality of 
markers. Preferably, data for the same set of maikers is available for each patient This 
set of markers may include all candidate markers, which may be suspected as being 
relevant to the detection of a particular disease or condition. Actual known relevance is 
not required. Enzbodiments of the methods and systems described herein may be used to 
detemune which of die candidate markers are most relevant to the diagnosis of &e disease 
or condition. 

[0085] The levels of each marker in the two sets of subjects may be distributed across 
a broad range, as illustrated in Figure 1 . Furfeer, although Figure 1 illustrates a 
distribution for the marker levels of the two sets, data for the two sets may simply be 
available as data points for each patient, as illustrated in Figure 2. No specific distributioa 
fit is required. 
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[0086] As noted above and as illustrated clearly in Figures 1 and 2, a marker often is 
inc^[»ble of effectively identifying a patient as either diseased or non-diseased. For 
example, if a patient is measured as having amaiker level that fells within the overlapping 
r^on, the results of the test may not be clinically relevant 

[0087] A cutoffmay be used to disting^iish between a positive and a negative test 
result for the detection of the disease or condition. Changing the cutoff trades off between 
the number of fake positives and the number of filse negatives resulting 6om the use of 
the single markier, or m the method described herein, the panel response. 

[0088] The effectiveness of a test having such an overlap is often expressed uang a 
ROC (Receiver Operating Characteristic) curve. OUksr measures, such as positive 
predictive value (PPV) and negative predictive value (NPV) may also be used as a 
measure of the esffectivencss of the test. ROC curves are well known to those skilled in 
tile art. For furtiier details, see Zweig, MH. & Campbell, C.C., Clin Chem 32, 561-577 
(1993) aad Hendrson. A.R., Ann. Clin. Biodhem 30, 521-539 (1993). 

[0089] Figure 3 illustrates an example of a ROC curve for the marker level 
distributions of Figure 1. The horizontal axis of the curve represents (1- specificity), 
which increases with the rate of false positives. The vertical axis of the curve represents 
sensitivity, which increases with the rate of true positives. ITius, for a particular cutoff 
selected, tiievahies of specificity and sensitivity may be determined. The area under the 
ROC curve is a measure of the utility of the measured marker level in tiie correct 
identification of one or more diseases or conditions. Thus, the area under the ROC carve 
can be used to determine tte effectiveness of the test. 

[0090] As discussed above, Ihe measurement of the level of a single maiker may have 
limited usefulness. The measurement of additional maikeas provides additional 
information, but die difficulty lies m properly combining the levels of two potentially 

unrelated measurements. 

[0091] In the methods and systems according to embodiments of the present invention, 
data relating to levels of various markers for the sets of diseased and non-diseased patients 
may be used to develop a panel of markers to provide a useful panel response. The data 
may be provided in a database such as Miciosofl Access, Oracle, other SQL databases or 
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siinply in a data file. The database or data file may contain, for example, a patieat 
iden^er such as a name or number, flie levels of the various markets present, and whether 
Ihe pati^ is diseased or non-diseased. Thus, a chart similar to Figure 2 may be generated 
for each marker of interest. Li practicei, the generation of Ihe diart is generally not 
required since the data may be directly acc^ble thioui^ Ihe database or the data file. 

[0092] In a preferred embodiment, one or more 'derived markers', which are a 
jftmction of one or more measured markers, may be incorporated into the set of markers 
being studied. For exaixple, derived markers may be related to the change m one or more 
measured marker values, or may be related to a ratio of two measured marker values, or 
may be other panel responses. In many diseases there will be a rapid change in a marker's 
vahie some time after an event. For exaoople, following an acute myocardial infarction, 
(AMI), myoglobin may rise rigidly and peak about 3 hours fiom the event. It may Ihen 
decay hsu± to its nominal value. Looking for changes in markers can be powerfiil 
diagnostic tool. Huis, &e change in myoglobin over a period of an hour, for example, 
may be used as a "marker" in the panel. 

[00931 ^ practice diagnosis of a disease state from multtple markers can be confusing. 
Often the individual maiker values may seem to contradict one another. In panels where 
the individual markeais are not very effective, it is extremely difScult to understand their 
meaning. In a preferred embodiment, a function that combines the marker values into a 
scalar value that increases with increasing likelihood of disease is defined. In this maimer, 
the informatiofn fixmi multiple markers may be presented in a useable form. This defined 
function is lefotxed to herein as the panel response (PR), and is a function of Ihe marker 
values (M^), written asPii = /(Jtfo-»)- "IT^fiwictional form of the pane^ 
be derived fi*om knowledge of the pathology of the disease, or may be part of the search 
space. Including the fimctional form as part of the search space may significantly increase 
the number of degrees of fi:eedom. In some cases the problem may become under 
specified. However, searching for the functional form may lead to further understanding of 
the pathology of the disease and the markers. The panel response may be scaled such that 
all values are between 0 and 1 . Because the effectiveness of the test may not d^end on a 
scaling of the panel response, scaUng may not influence the result of the method. However 
forcing ttie panel response to be a given scale may remove an unneeded redundancy, as 
panel response fimctions that differ only by a scaling factor may in fact represent the same 
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solution. The panel response may also be a general function of several parameters 
including tiie marker levels and o&er factors including^ fin: example, a patient's bistoty, 
age, race and gender of Hie patioat. 

10094] In a prefbrred embodiment, the panel response (PR) for each subject is 
expressed as: 

where i is the marker index, Wi is fhe weighting coeflBcient for the maricer i. Mi is the 
marker value for marker i, I is an indicator function for maiksr i, and X is the summation 
over all candidate madseas. The weighting Actors scale the indicator functions and may 
aUow for more important or q)ecific markers to have a greater imp^ 
response. The indicator fimction maps Hie marker value into a fimctianal fbtm appropriate 
to the marker's padbology. The indicator functions can be complex and sbould be chosen 
to match the marker. Hiis will be illustrated in the embodiments described below. The 
indicator fimction maybe a different fimctional form for each marker. In one example, the 
indicator function can m^ the marker value into a probability of the disease state. This 
mapping may not be a simple function of the marker value. In this exanq)le tiie said 
indicator fiom each marker can be summed to detemiine a relative index which is related 
to the probability of the patient being diseased. In a preferred embodiment the sum of all 
the weighting coefficients is constrained to a particular value, such as 1.0. As described 
above, this removes redundancy but doesn't change the objective function. In a preferred 
embodiment the indicator function is constrained to values between 0 and I . In a finUier 
preferred embodimoit, both of fhe above constraints are satisfied, thus, the panel response 
is also constrained to a value betweea 0 and 1 . 

{0095] In many disease states such as stroke, nonspeci&c markers associated with that 
state are elevated. But above a certain threshold, higher values of the marker may not 
relate to a higher probability of disease state. Below a certain threshold, lower marker 
v£dues may not relate to a lower probability of disease state. In this situation the indicator 
function may not increase linearly with the m^er value. A preferred embodiment is an 
indicator fimction that is a function that has a high and monotomc rate of change between 
the thresholds, and a snudl rate of change elsewhere. Examples of this type of function are 



20 



wo 2004/0580SS 



FCT/US2003/04i426 



the ramp, step, or sigmoid functions. One may associate the lower threshold with the start 
of an overlsq) tegion (or cutoff region), and the iq>per threshold with the end of the overlap 
region, as shown in Figure 4. Below the lower threshold the probability of disease is 
substantially 0, while above Ihe upper threshold flie probability of disease is 1 . Note that in 
the case wh^ l3ie indicator fvmction is a step function and the weighting value is 1 for 
each marker, tfien &e panel response is simply the number of markers above the cutoff. 
This case is identical to the exaniple used above where ono is searching for 1he best panel 
witii n of m markets above thdr cutoff. Allowing the indicatDr to vary continuously near 
ttie Hiieshold enables the pand response to be sensitive to a maiker just under the cutoff. 
This information is not lost as it is in Ifae n of m marker example or the step fonction 
example, ^ere the indicator value is not continuous. Anolher common approach of 
summing ever M* W forces the linear relation with the marker vahie» M. But as discussed 
above the most appropriate indicator function may not increase linearly with the marker 
value. In a further preferred embodiment the ramp function is used as an elevation 
indicator function. As illustrated in Figure 4, the indicator values betweeai the threshold 
regions may vaxy linearly jBrom a value of zero at one end to a value of one at the other 
end. In other embodiments, non-linear variations of the indicator value may be used. The 
ramp function has the advantage of simplicity, and may be good approximation to other 
fcmction in this class. With proper choices of parameters, the ramp function can be 
equivalent to the step function or can increase linearly with the tnaiker value. 

[0096] In some disease states, for example unstable angina, a ^lecific marker such as 
the cardiac troponins (including isofonns of cardiac troponin, comprising troponin I and T 
and complexes of troponin I, T and C) may be elevated above the normal population, but 
further elevation indicates an acute condition, in this case a myocardial infarction. 
Unstable angina is an ischemic condition that leads to minor necrosis of cardie tissue. 
During a myocardial infarction, there is major necrosis of cardiac tissue. Cardiac troponin, 
which is specific to cardiac necrosis, is elevated in hofti conditions, but the amount of 
elevation is related to the amount of necrosis. The best indicator function of cardiac 
troponin in diagnosing unstable angina may not be an elevation indicator function. In a 
preferred embodiment the indicator function may be a function that is peaked near the 
expected values of unstable angina, and decreases when the marker value is above or 
below the expected value. Exanqples of this type of fonction include a Gaussian, triaii^e, 
tapezxAd, or square fimction. These functions tend to localize fte maiker value of inteiBst 
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around a specific value. Anotfaer example of us© for such an indicator function is in cases 
where 8 pattern of markers values indicates a disease state. For example, a disease 
condition may be indicated when one or more markera are within a range of values. When 
desired, the use of this type of indicator may allow for recognition of patterns of maricer 
values. 

[0097] It is possible that one of the markers in the panel is spedjSc to the disease or 
condition being diagnosed. An example of sudi a marker is cardiac specific cardiac 
troponin when used in the diagnosis of acute myocaixlial infarction. The role of Tnl is 
described above. The panel response can be coded for markers that are specific, and the 
infbrmatiou may be used during the optimization of the panel response parametaB. 
Typically the specific cutoff of such maricare is known, so the specific cutoff values may 
not be included as a search parameter. When such a marker is present at above or below a 
certain flmeshold (e.g., the specific cutoff), Ihe panel response may be set to return a 
"positive" test result, regardless of the levels of non-apedfic matkas. When the specific 
cutoff is not satisfied, however, the level of the specific marker may nevertheless be used 
as a possible contributor to the panel response, along wifli the remaining markers on the 
panel. 

[00981 In an exanogjle whae the panel is being chosen based on n of m markers being 
elevated, the effectiveness of the panel is dependent on the choice of n. This extra 
dimensionality cm be eliminated by using an objective fimction. The reduction of 
dimensionality may shnplify the search process, and the objective function provides a 
scalar value that is optimized during the search process. The objective function should 
generally be indicative of the efiEectiveness of the panel, as may be expressed by, for 
example, overlap of the panel responses of the diseased set of subjects and the panel 
responses of the non-diseased set of subjects. In this maimer, the objective function may 
be optimized to maximize ttie effectiveness of the panel by, for example, minimizuig the 
overly. In a preferred embodiment, the ROC curve r^vresenting the panel responses of 
the sets of subjects may be used to define the objective function. A ROC curve witii a high 
value for the ROC curve area indicates a test with a good ability to discriminate between 
diseased and non-diseased. So, continuing with the n of m example above, ihere should 
exist a value of n which yields a dmically relevant test, but the value of n need not be 
detramined during the search process. The objective fimction is tiie scalar response that is 
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maximized by the search algoritibm. Hie objective fijnction can inclxide the correlation to a 
quantitative measure of disease, for example, the NIH score for stroke patients. Other 
measures of effectiveness may include, for example, odds ratios, a positive predictive 
value (PPV) and a n^adve predictive value (NPV) of fbs pand. The odds ratio, PPV and 
NPV are well knowu to those skilled in the art One skilled in the art will recognize there 
olher measures of the efiFectiveness of the test. See The Immunoassay Handbook, Second 
Edition, David Wild, 2001 for measures of efifectiveness. Many common measures of 
effectiveness require Ihe selection of a cutoff value. These functions may still be used, and 
Ihe cutoff value may also be included as a search parameter. In a prefisrred embodiment 
objective functions are chosen tibat do not require ihe selection of a cutoff value. The 
measure tibst is most aipropiiate for defining an effective test may vary. 

[0099] In a {neferred einbodimfint, the area uider the ROC curve lepresegating the 
panel responses of ttie sets of subjects may be used to define the objective fimctioo. Those 
skilled m the art will recognize that the area of the ROC curve is a measure of tiie 
effectiveness of Ihe test. An area of 1 corresponds to a perfect test, and an area of O.S 
corre^Kmds to a random test 

[0100] In another embodiment, the knee of the ROC curve is used for the objective 
function. The knee of the ROC curve is the point illustrated in Figure 6, and the value is 
represented as the product of the specificity and sensitivity at the knee, fii one embodiment 
the knee is found by maximizing the product of Specificity and Sensitivity. Higher knee 
values may mdicate squarer ROC curves. 

[01011 hi anothor embodiment the objective fimction is the specificity at a prescribed 
sensitivity. If one requires that a test have only a certain sensitivity (ability to detected 
diseased patients) tiien maximizing the specificity, which may reduce the number of false 
positives, may in^oye the clinical effectiveness of the test 

[0102] In anofliea: embodiment Ihe objective function is the sensitivity at a prescribed 
specificity. If one requires that a test have only a certain specificity (the number of false 
positives), then maximizing the sraisitivity, which may increase the abiUty to detect 
diseased patients, may improve the cUnical effectiveness of the test. 
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[01031 Th« objective fimction is a fimction of two or more of theae characteristics. In a 
preferred embodiment, the objective function is the pioduct of two or more characteristics 
of fihe ROC curve. An exarople of lliis is to use 1h« product of the ROC curve area, knee, 
sensitivity at a prescribed spedfidty, and specificity at a prescribed sensitivity. Any one 
characteristic alone may not result in a desired solution. By luhig the product of two or 
more of these, a more desirable solution maybe achieved. 

[0101] Variations in the values of markers over some time interval within a patient 
may be a powerful tool in the diagnosis of disease states or condition or the progression of 
disease states or conditions. As defined above die panel response can be thought of as a 
new derived marker, wbere the panel response value is thought of as the madcer vahie. 
Changes in the panel response value over some time interval within a patient may be a 
powediil tool in the diagnosis of disease states or conditions or the progression of disease 
states or conditions. The change in the panel response can also be used as a derived 
maiker. One can apply all of the ideas and methods discussed in this document to the case 
wh£!re a derived marker is a panel response or responses or the change in the panel 
response or responses. Calling the change in the panel response a derived marker may be 
equivalent to defining a new panel response tiiat is the change in the panel response over 
some time interval. The new panel response function is a functian of the marker values at 
two time points. All methods and ideas discussed in this documeait can apply to the new 
panel response. 

[0105] The two subject populations C'diseased" and "non-diseased") may contain 
subgroups. Thus, in accordance with the methods described herein, panel responses may 
be defined between such subgroups, considering one subgroup as a new "diseased" 
population and a second as a new "non-diseased" population. Accordingly, a panel 
r^onse combining the new panel responses for fliese new "diseased" and '"non-diseased" 
populations can be defined For exan:q>le, in a stroke data set, an origmal "non-diseased" 
population may contain both normals and stroke "mimics," Similarly, the "disesased" 
population may contain both ischemic and hemorrhagic stroke. Panel responses can be 
defined for normal vs. mimic, normal vs. ischemic, normal vs. hemorrhagic, mimic vs. 
ischemic, mimic vs. hemonhagic, and ischemic vs. hemorrhagic. Once these panel 
responses are defined, a panel response using these five panel responses for inputs either 



24 



wo 20047058055 



PCT/US2003/041426 



alone or in conjunction with Hie otbier maikas can be de£ba.ed to distingnish non-diseased 

fiom dise^ed. 

[0106] Searching for the best paael can be accomplished by trying all the different 
cambinations of parameters of the panel response function. But vrith panels of 40 markers, 
and just one degree of fteedom per marker, taking 10% steps in the parameter values will 
require lO'"* iterations. The age of the universe is estimated to be about 20 billion years or 
about 6.3X10^^ seconds. Clearly tins e^proach is not practical, and the problem requires 
the use of a search engine. Optimization algorithms are well-known to those skilled in the 
art and include several commonly available minimizing or maximizing ftmctions including 
the Simplex method and other constramed optimization techniques. It is understood by 
those skilled in the art tiliat some minimization fimctioDS are better tban others at searching 
for global nrniimunis, rather than local minimums. Many of these exist, and detailed 
descriptions can be found in the literature. For more infonnation on minimization and 
maxumzaHoin functions, reference may be made to Nnmeirical Recipes in C, The Art of 
ficf^tjfc nnyTipntinp. Second Mtiao, W. Press, et al, Cambridge UmveiBity Press, 1992, 
which is hereby incodpoiated by reference. Hie panel response and the objective function 
have he^ed enable Hie i^ of search routtnes. The objective function vahte is the response 
tbat the search routine will maximize, and the pazamdiers of the panel response fimction 
fbrm the n-dimensional space to be searched. While the objective fimction does not need 
to be continuous, Le. it may have discrete vahies, panel response functions that are 
continuous may reduce the granularily of the objective fimctioo. This may help the 
algorithm find better solutions. While many search routines will in fact look for minima, 
the problem may be inverted by minimizing (-l)*Objective Function. 

[01071 In a preferred embodiment the search engine uses the Downhill Simplex 
Method in MultkKmensions. This method is described in Numeripal Recipes m C. The Art 
of Scientific Computing. Second Edition, W. Press, et al., Cambridge University Press, 
1992. The simplex has iH-1 vertices, where n is the number of dimensions or degrees of 
fieedom. The routine 'walks* the siniplex along the n dimensional surface, moving one 
vertex at a time. The scale of the simplex can change so it can both quickly walk in 
downhill directions and crawl through tight crevices. The routine may not find a global 
minimum because it can become trapped in a local minimmn. The simplex will search all 
real space. The parameters of the panel response are often valid only within some range, 
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defining the bounds of tiie systeto. The sioiplex must be constrained to only search in Uiis 
space, and there must be no degeneracy introduced when api>ioacfajng such a constraint. 
One skilled in the art will recognize that tiKa» are many ways to address this consliaint In 
a preferred anbodimfflit a penalty is assessed when a vertex moves out of bounds. This 
penalty creates steep canyon walls around the bounds of the system, effectively 
constraining the amplex within the bounds of the system. 

101081 A weU-knownlimitatiom of search eaigmes is their tendency to 
minimum , typically not the global minimum. Several techniques are known to improve Ihe 
ability to seek out the global minimnm. In a piefsrred embodimeint, the technique of 
simulated annealing is used. This method is also described in ISTiimwri^jt^ Pg cioes in C. The 
Art of Scientific Cotnpntitig. Second Edition, W. Press, et al., Cambridge University Press, 
1992. Simulated annealing adds a random error to each decision of the search «igine. This 
random error gives the search engine the ability to move out of a shallow local minimum, 
so it can seek out a deeper solution. The random error is systematically reduced until a 
nunimum is found. The random earor is shnilar to the effect of temperature in annealing 
processes. The scale of the random error is said to be Ihe temperature. The annealing 
process may improve the chances of finding a global, rather than local, minimum, The 
annealing process may result in a more stable solution since the random variation may 
move the simplex out of a narrow, unstable region. The optimization process may be 
termioflted when the difference in the objective ftmction between two consecutive 
iterations is below a predetermined tiireshold, thereby indicating that the optunization 
algorithm has reached a region of a local minimum. The number of iterations may also be 
limited in the optimizatioa process. 

[0109] The selection of the initial conditions, for example the initial simplex value, 
may affect the optimization process. So, generally good selections of Uie initial parameters 
are sought. In the exair5)le of a search using a simplex, all vertices of the simplex must be 
initialized. If only one good vertex is defined, the other vertices can be assigned by 
flying a random deviation to each parameter. The scale of this random deviation sets the 
scale of the initial simplex. For example when elevation indicator functions are used, the 
location of the cutoff region may mitially be selected at any point But, selection near a 
suspected optimal location may &ciUtatB faster convogence of the optimizer. In a 
preferred method, the cutoff region is initially c«itered about the center of the overiap 
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region of the sets of patieaits. In one embodiment, the cutoff region may simply be a 
cutoff point In other embodiments, the cutoff region may have a length of greater than 
zero. In ihis regard, the cutoff region may be defined by a center value and a magnitude of 
lesDg&L la practice, the initial selection ofHae limits of the cutoff region may be 
determined according to a pie-selected percentile of each set of subjects. For example, a 
point above nvhich a pre-selected perceaitile of diseased patieaits ate measured may be used 
as tbe right (i^er) end of the cutoff region. In another embodiment the weighting &ctors 
may initially be all set to one. In a prefiBicred embodiment, the imtial weighting coefficient 
for each marker may be associated with the effectiveness of fiiat marker by itself. For 
example, a ROC curve may be generated for the single marker, and the area under the 
ROC curve may be used as the imtial weighting ooefficiait for that marker. This gives 
more weight to markers with better univariate utility. 

[QUO] Having seleotedoptunalparanieters&r the panel response function, the panel 
responses for each subject in each set of subjects, and the distribution of the panel 
responses for each set may now be analyzed. Figure 9 shows the ROC curves and area of 
several markers that have a poor diagnostic utility. The markers data are used to generate 
Figure S. When the poor markers are combined and the panel response determined, the 
results sbow that llie panel now has enhanced utUity. Figure S illustrates an exemplary 
distribution of the panel responses for diseased and non-diseased subjects. Based on these 
distributions, a ROC curve may be g«Dierate4 as illustrated in Figure 6. The ROC curve 
illuslrated in Figure 6 reflects optimized values for the weighting coefficients and the 
diresiholds for a ranq) mdicator function. 

[0111] Figure 7 illustrates an eatenq)lary progression of a ROC curve through a 
plurality of iterations of an optimization process in which the objective function is defined 
as the area under the ROC carve. As illustrated in Figure 7, as the number of iterations 
increases, the area undea: the curve may progressively increase. Thus, the optimization 
process may provide a panel response fonction for the markers. In this exanqjle, the 
indicator function is a ramp function. The optimization routine found values of llie 
weighting coefficients and high and low threshold values which are represented as a cutoff 
value and linear range. Table 1 illustrates a panel of 3 8 candidate markers with weighting 
coefficients and cutoff regions resulting fiwrn the optimization process. The 38 markers 
are listed generically as Analyte 1 tturough Analyte 38. The sense of each marker, as 
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described above, is also indicated in Table 1, with "liicr" represeating a positive sense and 
"Decsr" represoiting a negative sense. Hie cutoff location indicated in Table 1 refers to 
the marker level value aroimd viiiich the cutoff region is centerod, while the length of the 
cutoff indicates the range of marker level values covered by ttie cutoff region, hi tiiis 
manner, any number of markers may be used to develop a highly effective panel response 
function that can be used for the diagnosis of a disease or condition. 

[0112] The result of any given search is likely not to be the global minimum. It may be 
any tocal minimum that Ihe search engme settled in. In a product to be used &r clinic 
diagnosis, it is preferable to find a vcay stable rohition. loaocuracy associated with the 
measurement of the marker values should not significantly influence the efifectivesess of 
flie test. Also, the defiimg data may not be inclusive of all patients; it may be only a small 
sample, and the remainmg population may deviate from the defining population. The 
desired characteristics of the minimum may inchide a wide width and shallow side walls. 
In a three-dnnensional anrfogy, we would prefer a miniTmiTn Kke a crater as opposed to a 
mine shaft. One method to seek out fliese types of solutions is to search multiple times. If a 
statistically significant number of optimizations is performed, then the better solutions will 
be the largest group of similar results. This is because, uaing the example above, it may be 
more Hkely to find tiie crater tiian the mme shaft 

[0113] As discussed above, not eveary TniTiim mfn fbund may be desirable to use. 
Generally stable paramet^ are desked, meaning tiiat variations in Ihe marker values or 
parameters do not adversely impact the effectiveness of the test. The width and deptii of 
the minimum may provide an indication of the stabiUty of the solution. Li addition, one 
could use bootstrapping methods known in the art to determine the stabiHty of a particular 
solution. In such methods, a "training" data set is used to arrive at an initial solution. That 
solution is then, validated by applying the solution to a second "validation" data set that is 
independent of tiie traming data set 

[0114] There are several additional exanqples of methods that may quantify the quality 
of a set of parameters, e.g., whsa a validation data set is unavailable. A first example is to 
vary the maifcer values by some random parentage. By doing this one can simulate all the 
variations expected due to assay imprecision, biological variations, and any otiier source 
of unceortainty. For example, variatioojs in marker values may relate to the relative 
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imprecision of the test tiiat was used to gaiCTate ttie data. One skilled in the are will 
recognize timt there are liniits to the analytical precision of a test. For example, in the 
immunoassay art, it is common to encounter 5-20% coefficients of variations of the tests. 
Thetefore, when considering the imprecision of the testing methodology, the paiameters 

stable relative to the imprecision of the m^Jiodology. The randomized data set can 
be reanalyzed to generate Hie new panel response ROC curve and objective function value. 
An acceptable deterioration may indicate the parameters give a solution that is stable to 
variations in maiicer values and may also verify that the solution does not siniply fit liie 
noise in the data. 

[01151 A. second example would be to vary one or more of the parameters in the panel 
definition by some amount. The change in the objective function value may be a measure 
of the quality of the solution. Each parameter could be varied independently to detemine 
the stability of each parameter. Or all parameters may be varied randomly to sample the 
space surrounding the solution. The magnitude of the variation can be constrained to 
ssaapls a certain volume or shell surrounding the solution. Thus the deterioration of the 
objeiCtive function as a function of displacement from Ifae sohition maybe detmnined. 

[0116] In a third example, a seed simplex is generated with a given length scale about 
liie known mini mum. The length scale of the seed simplex can be systematically increased 
until re-optimizations lead to a dif^ent minimum, i.e. the solution is no longer recovered. 
The length scale, which results in finding new mlDimums, may be related to the width of 
the minimntn. In a fourth example, using the final simplex of the optimization, the 
tenqjerature can be systematically increased until re-optimizations lead to a different 
minimum, i.e. the solution is no longpr recovered The temperature, which results in 
finding new mininmms, maybe a measure of the depth of the solution. In a fiftb example, 
most common solutions fium the multitude of optimizations, may represent the most 
stable solution. The common solutions can be grouped based on their similarity. 
Correlation techniques and clustering techniques can be used to gcaap the solutions, and 
are well known to one skilled in tibe art. Solutions can also be grouped by fitKling ones 
that, using a three dimensional analogy, are on the same valley floor. These can be 
grouped by calculating the reduction in objective function while movmg ftom one solution 
to another. In the analogy above this detennines how high of a hill lies between the two 
solutions. Fixnn the teaching above, it is now clear that other appioadies exist for 
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quanttfying flie quality of a set of parameters, and the examples above are not intended to 
limit the invention. 

101171 The use of the term "non-diseased" does not mean that the particular subject is 
disease-free, only that the subject is free from the one or more diseases or conditions being 
evaluated In practice, a pre-filtering of subjects may be performed on the basis of any 
particular characteristic of the subjects, including the existence of other diseases. For 
esxaniple, the method and systems described may be applied to first divide a groi^ of 
subjects into "diseased" and "non-diseased" for Disease A, and thai divide liie group into 
"diseased" and "iion-diseased" for Disease B. A panel of madcets for each disease may 
then be determined. In another embodnnent, the same panel of markers may be used for 
both diseases with a diffaent set of parameters, such as weighting coefficients, for each 
disease. In another embodiment subjects with disease A can be defined as non-diseased, 
and subjects with disease B can be defined as diseased. In this embodiment the described 
techniques can be employed to determine a panel that differentiates between disease A 
andB. 

[0118] The search routine will optimize the objective function or fonctions selected on 
the specified data set But often times it is important to constrain or optimize a second 
group of data simultaneously. This is accomplished by pre-filtering the source data to get 
the two or more groups of data of mterest Different objective functions can be selected for 
each groiq) of data, and the search engme can find the mininwrn of the product of 
objective iiinctiQinB. The objective fimction of one of the groui^ of data can also be 
constrained to be at least some value. When the objective function is greater than or equal 
to this constraint, tiie value returned to the search engine is the constraint value. When the 
objective fimction is below the COTistraint value the objective function value is returned. 
The search routine will look for solutions that satisfy the constraint condition, but the best 
solution may fall outside the constraint condition. The iterations of the optimization 
algorithm generally vary the indq)endent parameters to satisfy the constraints while 
maximizing the objective function. An example of this usage is stroke data that contains 
norm health donors and stroke mimics. We would like to find a panel response function 
that will distinguish shoke from stoke mimics, but tiiat will also have a low false positive 
rate for normal heaMiy donors (NHD). Since the number in each sample set is not equal, 
siniply combining the data and analyzmg will not give a satisfectory result Results will be 
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skewed to the data set with larger n, in our case NHD. However, if the objective functions 
of the two groups of data are individually calculated and combined, tihen the groi5)s of 
data aie given equal weight. In another example we want to ensure that patients presenting 
soon after flie onset of symptoms will be propedy diagnosed, but we still want to eosure 
ikat patients preseniing at longer times are also properly diagnosed. Again, the population 
numbers will be different. So, to give equal weigbtuig. fhey need to be simultaneously 
analyzed as two gxnsps of daita. Other constraints may include limitations on one group of 
samples mlule optimiziiig for an objective function for a second group. For exan^le, a 
panel may be optimized for one disease wbile the same panel may be constrained to 
provide at least an acceptable minimum value for the area under a ROC curve for a second 
disease. 

[0119] Within the teachings of this document we have used for simplicity maikers lhat 
are elevated hi patients with the disease or positive sense marlcets. However this is not 
always the case, and often, particularly with pooar univariate maikers, it is not dear fram 
univariate analysis whether the marker when used in conjunction with the other maikers in 
the panel, is best utilized in a positive or negative sense. If the sense of a marker is 
inverted, then it is stisi^tfoxward to invert tiie indicator function for that marker. If the 
sense is not known, thai the search engine may include this as a degree of fieedom. For 
cxmaple, 'm one embodiment, the sense may be a truly separate independent variable, 
which may be flipped between positive and negative by the optimization process. For 
optimal perfonnattce, the sense should map smoothly firom hnproper to proper, and there 
sihould be pressure (or a gradient) that allows the search engine to move toward the proper 
sense. In a preferred embodiment the sense is switched by allowing the weighting 
coefficient of the analyte to go negative. If the wrong sense is selected, the weighting 
coefficient will be driven towards zero since inclusion of the marker in the panel response 
negatively in^acts the objective function. The search engme will be able to drive the 
wdghtmg coefficient across zero to the proper sense. In tiiis example, the negative weight 
is just a flag to invert the sense. Hie absolute value of Ihe weigiht is used in the panel 
response function. But tiiis allows for a continuous function moving firom positive to 
negative sense. 

[0120] In order to determine lije best panel, which for practical reasons may often 
mean 10 or less maikers, one must find a way to systematically remove markers that do 
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not significantly contribute to iJie overall result Again, there are sevctal methods that may 
be ^plied to arrive at a useful set of maikexs. For exabq>le, markers may be initially 
selected by performing univariate statistics to detexmine if a marker provides a meaningful 
distmciion between the "disease" and *^on-disease" groups. While the distinctiDn m such 
univariate methods may be weak, any correlation to the disease or non-disease state may 
be enough to indicate that a madcsr should be considered for further analysis. 

[0121] One may also calculate the contribution &om each nmker to Hie marker panel 
result itself. A method to accomplish this is to remove an analyte fiom the panel, and 
recalculate the objective function. This can be achieved by re-optimizing the parameters in 
the absence of the analyte and detennining the best objective fbnction. One can also 
remove the analyte and recalculate the objective function without re-optimizing. While not 
as precise^ it can offer significant savings computationally. In eitlier case, care must be 
taken to ensure that the data set is static -vAim the analyte is removed. The change in the 
objective function is related to the contributioii of the madcer. This method for identifying 
the relative importance of each marker is ilhistrated in Figure 8. The resulting changes in 
the objective fimction are noted for each marker and plotted, as shown in Figure 8. Figure 
8 illustrates the effect each mstkei has on the various features of the ROC curve 
corresponding to the panel responses for the two sets of subjects. Ite various ROC-curve 
features noted in Figure 8 include the area under the ROC curve, the location of the knee 
of the ROC curve, the sensitivity at a predetermined specificity, and the specificity at a 
predetennined sensitivity. The markers may then be arranged in order of decreasing 
contribution, as illustrated in Figure 8. The vertical axis in Figure 8 indicates tiie relative 
change in the values of tiie various ROC-curve features. In embodimeiits where a 
weighting coefScioit is applied to each analyte. the weight for the analyte can be set to 
zero to remove the analyte's contribution In tiie panel result. The change in objective 
function can then be determined. In embodiments whax a weighting coefQcient is applied 
to each analyte, one can not shnply use the weights as the contribution. An example of 
why this does not give die proper result is the case where a marker has zero impact on the 
test. In this case, the wei^ it is given by the search program can be any value, so it is 
po^ble that its weight will be the highest. 

[0122] In order to develop lower-cost panels, which require the measuremcsnt of fewer 
marker levek, certain markers may be eliminated fiom die panel. In this regard, the 
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effective contributioEL of each maiker in the panel may be determined to identify the 
lelative importance of the madkers. Once the lelative contributions are calculated thai one 
can rank Ihem fiom largest to lowest. The markers with the largest changes in objective 
fimctiQn may be the ones with most coolxibution. The ones with the le^t change in 
objective function may be the ones with the least contribution. However, for example^ if 
two markers are perfectly correlated, then the combined oonbibution fbom boHi may be 
equivalent to Ifae contribution of just one if the second one is removed. The partitioning of 
the contributions is not necessarily equal So an important marker may not have a high 
contribution. This problem can be avoided by first looking at the correlation, or 
'^interaction," between markers, or by removing only one marker or more witii the bwest 
contribution. Methods for calculating interaction terms are well known to those of skill in 
tiie art. 

[0123] From the discussion above, it is noted tiiat it may not be prudent to just select 
the top 3 maikers fiom a panel of 40. Depending om the nuniber of target madcers being 
seaidhed and the size of the target panel, one may want to eliminate only the marker with 
the lowest contribution or the lowest markers, and rq>6at tiie piocess until the target panel 
size is reached. With properly defined panel responses, markers of no importance may not 
adversely impact the objective function. This is because a) the search routine may chose 
parameters such that the marker is not used, and h) in general a random marker will not 
change the objective fimctioiL So, starting with a large panel and reducing it to the desired 
size will lead to the optimum panel. But the objective function may degrade as markers are 
eliminated. One may have to trade off panel effectiveness with ttie number of markers. For 
example, in order to obtain a panel of ten maik^s, the ten highest-rated markers, i.e. those 
on the left side in Figure 8, may be selected. For example, Analytes 38, 1, 16, 33, 27, 12 
and 8 may be selected in a final panel of markers. In a preferred embodiment, only a few 
of the markers on the right side may be eliminated, and the remaining markers in the panel 
maybe optimized. For example, Analytes 31, 24, 25, 4 and 10 maybe eliminated in a first 
round, and I3ie optimization and ranking procedures may be repealed with the remaining 
33 markers. This results in a chart similar to that shown in Figure 8, but witii fiswer 
markers. This process may be repeated until a desired number of markers remains in the 
panel. 
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[012^ It is possible fbat one of th^ markers in fhe panel is specific to the disease or 
condition being diagnosed. An example of such a marker is cardiac specific Tnl when 
used in tihie diagnosis of acute myocardial infarction. The role of Tnl is described above. 
The panel response can be coded for markers that are specific, and the information is used 
during liie optimization of the panel response parameters. Typically the cutoff of such 
markers is known, so tiie cutoJff values may not be included as a search parameter. When 
sudi a marker is present at abovo or below a certain threshold, the panel re^onse may be 
set to return a '^positive" test result, regardless of the levels of non-specific maricers. 
When the threshold is not satisfied, however, the level of the specific mark^ may 
nevertheless be used as possible contributor to flie objective fiinction, along with the 
remaining markers on Ihe panel. 

[0125] In a prefmed onbodiment the panel will include markets derived fix)m the rate 
of change of markers measured by the panel. In a fortiiiOTprefeired embodiment the panel 
will have two panel response functions, one that utilizes tiie derived markere when present, 
and when not presait one that does not utilize the derived markers. The two panel 
response functions may use different parameters. These parameters may be obtained by 
optimizing the data with and wittiout utilizing the derived marker or markers. For 
example, a patient may be measured when first arriving at the hospital for a particular set 
of markers. Smce there is only one sample time for the patient a panel response fimctioai 
which does not include marker changes is used. The patient would be diagnosed as 
diseased <x non-diseased based on the results of the test. The same patient may be 
measured again an hour later. Now th^ axe two pomts, and so a second panel response 
function which utilizes marker changes is used. Hie use of this response function is 
important vfbm a marker or panel of markers of disease indicates non-disease, but the 
change (usually increase) in fhe value of one or more mariceis represents Ihe onset of 
disease. 

I0126I It is possible for a panel of markers to contain enough information to diagnose 
a multitude of conditions. In the simplest case, the markers used m the diagnosis of 
condition A are different fiom the markers used in the diagnosis of condition B, but the 
panel contains Ihe union of fhe markers for A and B. In a prefeired embodiment, the 
markers used in the diagnosis of condition A contains at least one of the markers used in 
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the dis^osis of condition B. lii a future prefiEiTed embodiment there is a high degree of 
overlE^ in the markers used to dia^ose a inuttitude of conditions. 

[0127] The method described above may be implemeaited in a variety of manners. Lia 
preferred embodiment, the method is implsnoited as a program product, such as a 
software package. The program product may be implemented on a computer, such as a 
personal computer, a main&ame or a handheld device. It will be apparent to those skilled 
in the art fliat the program product may be implemented on a device m any number of 
ways including software, firmware, etc. In one embodiment, the program product is 
implemented on a meter which may be capable of directly measuring leveis of one or more 
markers. For exaitq)le, the program product may be implemented on a fiuoiometer or a 
reflectoraeter. Such devices are well known to those skilled in the art. 

10128] In a most preferred mode, patient types, disease types, and time frames are 
selected to provide two data sets, "diseased" and '*non-diseased," which have the 
characteristics to be evaluated. Multiple groups of data can be selected, each set consisting 
of a set of diseased and as a set of non-diseased samples. The values for any derived 
marker values of interest are calculated for each record in the selected groups of data. This 
may include calculating the change in marker value or panel response value fix>m tbo 
initial vahie. Based on llie disease and marker pathology, a functional type for the 
indicator function is chosen for each maiker to be included in the panel The teadhings in 
tiiis document should enable one skilled in the art of the disease and maiker to make fins 
apptopiiaite dioice. if the disease or maiker pathology is not sufficiently understood to 
choose a functional form th^ the functional fbnn can become part of the search. Once liie 
indicator functions have been defutied, Iheai the initial parameters ace chosen &am the 
univariate marker analysis. These initial parameters define one vertex of the initial 
sing)lGX. The number of vertexes constituting the simplex is the number of search 
parameters in the panel respoi^e phis one. Each remainmg vertex is populated by varying 
each parameter by aiandom amount The scale of this random amount can be fixed to be a 
percentage of Hie parameter value. Hub spreads the simplex out around the initial point, 
and gives the simplex a size scale. The objective function for each groiq) of data is defined 
by selecting any combmation of the ROC curve area, the ROC knee, the ROC sensitivity, 
and the ROC specificity, but typically all four are selected. The objective fimction for each 
group of data can be chosen to be optimized or to Tnainffliti s, minimum target value. Thus 
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the optiioization of one groop of data can be constrained such that a second groiq) of data 
has at least a miTinmim objective function value. The parameter ate then optmdzed to 
maximize tfie chosen objective function utilizing Ihe downhill simplex msthod with 
simulated amiealing. At the end of the optimization Ihe relative contribution for each 
maiker is calculated by setting the weight of lhat marker to zero and recalculating the 
panel ROC curve. When the analyte is so removed from the panel response, the new ROC 
curve is calculated with fbs identical data and no other parameters in the panel response 
are changed. The process of optimizing and calculating marker contributions is repeated n 
(~100) times. After n optibraizations, the averse contribution of each marker over the n 
optimiza:tioiis is calculated, and the markers are ranked based on its average contribution. 
The poorest markers, typically the poorest half or less, are removed from the panel and the 
entire process is repeated as many times as required to reduce the panel to the desired size. 

[0129] Candidate solutions are then tested for stability and robustness by varying the 
input data and the paiametera, and finding solutions that are groined together. The best 
solution may be derived fix)m an average of one or more good solutions. Using optimal 
analytes and parameters for ttie panel response fimction found via the search method 
described above, tiie ROC curve of the panel response from clinical data is calculated. 
Based upon the panel response ROC curve an appropriate cutoff is chosen. The choice 
maybe influenced by factors such as cUnicd factors, treatment methods, and cost 
considerations, which one skilled m the art will recognize. The panel response is 
calculated from the measured marker values of the patient for whom it is desired to 
detenmne the presence or absence of the target disease. Using the chosen cutoff assign a 
diagnosis for the patient. 

[0130] Using optimal analytes and parameters for the panel response function found 
via ttie search method described above, for panel response functions which include and 
exclude markers derived from the change in a measured marker, the ROC curve of the 
panel response from clinical data is calculated. Based upon the panel response ROC 
curves ^piopriate cutoffe are choseaa for each. The choice may be influenced by factors 
such as cUnical factors, treatment methods, and cost considerations, which one skilled in 
the art will recognize. Upon measuremrait of the initial sample, the panel response is 
calculated from the measured marker values of the patient for whom it is desired to 
detemiine itss presence or absence of the tai^ disease. Using tiie chosen cutoff assign a 
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diagnosis for the patient. A second or more measurement may be required to further 
clarify the diagnosis. At the appropriate time interval, draw more sample from the patient 
and measure the marker values. Using the panel re^onse function that includes derived 
markers, calculate the panel response value and determine a diagnosis by comparing the 
panel response value to the chosen cutoff value. The panel response of tiie first 
measurement can also be compared to panel responses determined fiom subsequent 
measurements. One skilled in the art will recognize that serial blood draws can yield 
critical information of the presence and progression of diseases, particularly acute 
diseases. If more measurements are required for propffl: patient treatmsat, continue taking 
samples at the desired intervals. 

[0131] Measures of test accuracy may be obtained as described in Fischer et al.. 
Intensive Care Med. 29: 1043-51, 2003; ZhaaetaL, Statistical Methods in Diagnostic 
Medicine, John Wiley & Sons, 2002; and Motnlsiky, Intuitive Biostatistics, Oxford 
University Press, 1995; and other publications well known to those of skill in the art, and 
used to detennine the effectiveaiess of a given marker or panel of markers. These measures 
inchide sensitivity and specificity, predictive values, likelihood ratios, diagnostic odds 
ratios, hazard ratios, and ROC curve areas. As discussed above, suitable tests may exhibit 
one or more of the following results on Ihese various measures: 

[0132] A ROC curve area of greater than about 0.5, more preferably greater than about 
0.7, still more preferably greater than about 0.8, even more preferably greater than about 
0.85, and most preferably greater than about 0.9; 

a positive or negative likelihood ratio of at least about 1 , 1 or more or about 0.9 1 or less, 
more preferably at least about 1.25 or more or about 0.8 or less, still more preferably at 
least about 1 .5 or more or about 0.67 or less, even more preferably at least about 2 or more 
or about 0.5 or less, and most preferably at least about 2.5 or more or about 0.4 or less; 

an odds ratio of at least about 2 or more or about 0.5 or less, more preferably at least about 
3 or more or about 0.33 or less, still more preferably at least about 4 or more or about 0.25 
or less, even more preferably at least about 5 or more or about 0.2 or less, and most 
preferably at least about 10 or more or about 0.1 or less; and/or 
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a hazard ratio of at least about 1.1 or more or abotit 0.91 or less, more preferably at least 
about 1.25 or more or about 0.8 or less, still more preferably at least about 1.5 or mote or 
about 0.67 or less, even more preferaibly at least about 2 or more or about 0.5 or less, and 
most preferably at least about 2.5 or more or about 0.4 or less. 

[0133] Measures of diagoostic accuracy such as those discussed above ate often 
reported together with con£ideaice intervals or p values. These may be calculated by 
methods well known in the ait See, e.g.. Dowdy and Wearden, Statistics for Research, 
John Wiley & Sons, New Yori£, 1983. Prefcaied confidence intervals of the invention arc 
90%. 95%, 97.5%, 98%, 99%, 99.5%, 99.9% and 99.99%, while prefeied p values are 
0.1. 0.05, 0.025, 0.02, 0.01. 0.005, 0.001, and 0.0001. 

[0135] Patients presenting for medical treatment often exhibit one or a few primary 
observable changes in bodily characteristics or functions that are indicative of disease. 
Often, these "syoqitoms" are nonspecific, in that a number of potential diseases can 
present the same observable symptom or symptoms. A typical list of nonspecific 
symptoms might include one or more of the following: shortness of breath (or dyspnea), 
chest pain, fever, di2ziness, and headache. These syn^ms can be quite common, and the 
number of diseases that must be considered by the clinician can be astoundingly broad. 

[0136] Taking shortness of breath (refored to clinically as "dyspnea") as an example, 
this symptom considered in isolation maybe indicative of conditions as diverse as asthma, 
chronic obstructive puhnonary disease ("COPD"), tracheal stenosis, obstructive 
endobroncheal tumor, pulmonary fibrosis, pneumoconiosis, lymphangitic carcinomatosis, 
kyphoscoliosis, pleural effusion, amyotrophic lateral sclerosis, congestive heart failure, 
coronary artery disease, myocardial in&rction, cardiomyopattiy, valvular dysfunction, left 
ventricle hypertrophy, pericarditis, arrhythmia, pulmonary embolism, metabolic acidosis, 
chronic brondutis, pneumonia, anxiety, sepsis, aneurismic dissection, etc. See, e.g., 
Kelley's Textbook of Internal Medicine, 4* Ed., Lippincott WilKams & WiDdns, 
Philadelphia, PA, 2000, pp. 2349-2354, "Approach to the Patient With Dyspnea"; Muhx>w 
etal, J. Gen. Int. Med. 8: 383-92 (1993). 
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[0137] Similarly, chest pain, whai considered in isolation, may be indicative of stable 
angina, unstable angina, myocardial infarction, musculoskeletal injury, cholecystitis, 
gastroesophageal reflux, pulmonary emboMsm, pericarditis, aortic dissection, pneumonia, 
anxiety, etc. Mmeoyec, the classificaticm of chest pain as stable or unstable angina (or 
even mild myocardial infarction) in cases othei than definitive myocardial infarction is 
completely subjective. The diagnosis, and in this case the distinction, is made not by 
angiography, yUAch may quantify the degree of arterial occlusion, but rather by a 
physician's interpretation of clinical symptoms, 

10138] Differential diagnosis refera to methods for diagnosing the particular disease(s) 
underlying the symptoms in a particular subject, based on a comparison of the 
charactmstic features observable fiom the subject to the characteristic features of those 
potential diseases. Depending on the breadtii of diseases that must be considerod in the 
diffaxntial diagnosis, tiie types and number of tests that must be ordered by a chnician can 
be quite large. In the case of dyspnea for esisample, the cliniciaa may order tests fitm a 
group that includes radiognphy, electrocardiogram, exercise treadmill testing, blood 
chemistry analj^s, echocardiography, bronchoprovocation testing, spirome^, pulse 
oximetty, esophageal pH monitoiing, laryngoscopy, computed tomography, histology, 
cytology, magnetic resonance imagmg, etc. See. e.g., Morgan and Hodge, Am. Fam, 
Physician 57: 711-16 (1998). The cUnician must theai mtegrate information obtained from 
a battcary of tests, leading to a clinical diagnosis tiiat most closely represents the range of 
syn^toms and/or diagnostic test results obtained for tiie subject 

[0139] A first step in the identificatioa of suitable markers for symptom-bases 
differential diagnosis requires a considaration of the possible diagnoses that may be 
causative of tiie non-specific sycaptam observed. Takmg dyspnea as an example, the 
potential causes are myriad. The followutg discussion considers tisree potential diagnoses: 
congestive heart failure, pulmonaty embolism, and myocardial infaiction; and three 
potential markers for inclusion in a differential diagnosis panel for these potential 
diagnoses: BHP, D-dimer, and cardiac ttopooin. 

[0140] BNF 

[0141] B-type natriuretic peptide (BMP), also called brain-type natriuretic peptide is a 
32 amino acid, 4 kDa peptide tiiat is uivolved in die natriuresis system to regulate blood 
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pressure and fluid balance. Bonow, R,0., Circulation 93:1946-1950 (1996), The precursor 
to BMP is synfliesized as a 108-anmio acid molecule, refetred to as **pre pro BMP," that is 
proteolytically processed into a 76-aminD acid N-tenninal peptide (amino adds 1-76), 
refiEwedtoas "T^iT pro Bl^" and the 32-amino acM mature homonci referred to as B]^ 
or BNP 32 (amino adds 77-108). It has been suggests that each of these species - NT 
pro-BNP, BNP-32, and the pre pro BNP - can circulate in human plasma. Tateyama et al, 
Biochem. Biophys. Res. Common. 185: 760-7 (1992); Hunt et al., Biockem. Biophys. Res. 
Conamm. 214: 1175-83 (1995). The 2 forms, pre pro BNP and NT pro BNP, and peptides 
vMch are daived fiom BNP, pre pro BNP and NT pro BNP and which are present in the 
blood as a result of proteolyses of BNP, NT pro BNP and pre pro BNP, are collectively 
described as markers related to or associated with BNP, 

10142] The tain 'BISEP'' as used herdn refers to the inatur©32-aminoad^ 
molecule itself As the skilled artisan will recognize, however, because of its relationship 
to BNP, the concentration of NT pro-BNP molecule can also provide diagnostic or 
prognostic information in patients. The phrase "marker related to BNP or BNP related 
peptide" refers to any polypeptide that originates firom the pre pro-BNP molecule, other 
tiian the 32-ammo acid BNP molecule itself. Proteolytic degradation of BNP and of 
pq)tides related to BNP have also been described m the Hterature and these proteolytic 
fragments are also enconipassed it the teon "BNP related peptides." 

[01431 Blyff and BNP-related peptides are pradommantly found in tiie secretory 
grannies of the cardiac ventricles, and are released from the heart m response to both 
ventricular volume expansion and pressure overload. Wilkins, M. et al. Lancet 349: 1307- 
10 (1 997). Elevations of BNP are associated with raised atrial and puhnonary wedge 
pressures, reduced ventricular systolic and diastohc function, left ventricular hypertrophy, 
and myocardial infarction. Sagnella, G.A., Clinical Science 95 : 5 1 9-29 (1 998). 
Furthermore, there are numerous reports of elevated BNP concentration associated with 
congestive heart failure and renal faihtte. Ihvs, BNP levels in a patient may be mdicative 
of several possible underlying causes of dyspnea. 

[0144] D-dimer 



[0145] B-dmer is a crosslinked jShrin degradation product with an approximate 
molecular mass of 200 kDa. The normal plasma concentration of D-duner is < 150 ng/ml 
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(750 pM). The plasma concentration of D-dimer is elevated in patients with acute 
myocardial infarction and unstable angitia, but not stable angina. Hoffineister, HJM. et al. 
Circulation 91: 2520-27 (1995); Bayes-Genis, A. et cd., Thromb, Haemost. 81: 865-68 
(1999); Gurfmkel, E. etal., Br. Heart J. 71: 151-55 (1994); Kruskal, JJB. et al., N. Engl. J. 
Med. 317: 1361-65 (1987); Tanaka, M. and Suzuki, A, Tkromb. Res. 76: 289-98 (1994). 

[0146] The plasma concentration of D-dimer also will be elevated during any 
condition associated with coagulation and fibrinolysis activation, including stroke, 
surgery, atherosclerosis, trauma, and thrombotic thrombocytopenic purpura. D-dimer is 
released into the bloodstream immediately following proteolytic clot dissolution by 
plasrain. The plasma concentration of D-dimer can exceed 2 ixg/ml in patients with 
unstable angina. Gurfinkel, E. et al., Br. Heart J. 71: 151-55 (1994). Plasma D-dimer is a 
specific marker of fibrinolysis and indicates the presence of a pioHirombotic state 
associated with acute myocardial infarction and unstable angina. The plasms conceaatratioaa 
of D-dimer is also nearly always elevated in patients with acute pulmonary embolism; 
thus, normal levels of D-dimer may allow the exclusion of pulmonary embolism. 
Egetraayer et al., Thorax 53: 830-34 (1998). 

[01471 Cardiac Troponin 

[0148] Troponin I (M) is a 25 kDa inhibitory element of the tit)ponin complex, found 
in muscle tissue, Tnl binds to actin in the absence of Ca^^ inhibiting the ATPase activity 
of actomyosin. A Tnl isoform that is foimd in cardiac tissue (cTnl) is 40% divergent from 
skeletal muscle Tnl, allowing both isoforms to be iinmxmologically distinguished. The 
normal plasma concentration of cTnl is < 0.1 ng/ml (4 pM). cTnl is released into the 
bloodstrefflu following cardiac cell death; thus, the plasma cTnl concentration is elevated 
in patients with acute myocardial infarction. Investigations into changes in the plasma 
cTnl concentration in patients with unstable angma have yielded mixed results, but cTnl is 
not elevated in the plasma of individuals with stable angina, Bemamer, H. c/., Am. J. 
Cardiol 82: 845-50 (1998); Bertinchant, JJ. et al., Oin. Biochem. 29: 587-94 (1996); 
Tanasijevic, M.J. et al, Clin. Cardiol. 22: 13-16 (1999); Muaao, P, etal, J. Ital. Cardiol. 
26: 1013-23 (1996); Holvoet, P. et al., JAMA 281 : 1718-21 (1999); Hohroet, P. et al.. 
Circulation 98: 1487-94 (1998). 
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[0149] The plasma concentration *)f cTnl in patients with acute myocardial inferction 
is significantly elevated 4-6 hours after onset, peaks betweei 12-16 hours, and can remain 
elevated for one week. Tho release kineticB of cTnl associated with unstable angina may 
be similar. The measurement of specific forms of cardiac troponin, inchiding fiee cardiac 
troponin I and complexes of cardiac troponin I witii troponin C and/OT T may provide the 
user with fee ability to identify various stages of ACS. Pree and complexed cardiac 
troponin T may be used in a manner analogous to that described for cardiac troponin I. 
Cardiac troponin T complex may be usefiil ei&er alone or when expressed as a ratio with 
total cardiac troponin I to provide information related to the presence of progressmg 
myocardial damage. Ongoing ischemia may result in the release of the cardiac troponin 
TIC complex, indicating that higher ratios of cardiac tropomn TICitotal cardiac troponin I 
may be indicative of continual damage caused by unresolved ischemia. See, U.S. Patent 
Nos. 6,147,688, 6,156,521, 5,947,124, and 5,795,725. 

[0150] Based on the foregoing discussion, the skilled artisan will recognize that, for 
example, increased BMP is indicative of congestive heart failure, but may also be 
indicative of other cardiac-related conditions such as myocardisd infarction. Thus, the 
inclusion of a marker related to myocardial injury such as cardiac trq)omn I and/or cardiac 
tropomn T can permit further discrimination of the disease und^lying the observed 
dyspnea and the increased BNP level. In this case, an uwaeased level of cardiac troponin 
may be used to rule in myocardial infarction. 

[0151] Similarly, BMP may also be indicative of puhnonary embolism. The inclusion 
of a matka: related to coa^ation and hemostasis such as D-dimer can permit further 
discrimination of the disease underiying the observed dyspnea and the increased BNP 
level. In this case, a normal level of D-dimer may be used to rule out pulmonary 

embohsm. 

[0152] The skilled artisan will readily acknowledge lhat other markers may be 
sut^tituted in or added to tiiis marker panel to further discriminate the causes of dj^pnea. 

Suitable markers are described in co-pending PCT Application No. ^ filed 

December 23, 2003 (Atty. Docket No. 071 949-5603), which is hereby mcoiporated by 
reference in its entirety. Preferred panels for the diagnosis of a cause of dyspnea comprise 
a phuality of markers independently selected ftom the groiq) consisting of specific 
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markers of cardiac injury, specific markera of neural tissue injury, non-specijEic markeis of 
tissue injury, mark^s related to blood pressure regulation, markers related to 
infUmunatioD, markers related to coagulation and hemostasis, markers related to 
pulmonary injury, and majjcers related to ^optosis. £x.eiiiq)lary markers in each of these 
groiqts are described heareinsfier. Preferably, such a panel comprises maikeaxl fiom two, 
three, four, jBve, or more difGsent members of tiiis group. Thus, particularly preferred 
panels for the diagnosis of a cause of dyspnea comprise one or more speciSc maikers of 
cardiac injury and one or more maikais related to blood pressure regulation; one or more 
spedfLc maikcrs of cardiac ingury and one or more markers related to coagulation and 
hemostasis] one or more matkeais related to blood pressure regulation and one or mcxe 
maikers related to coagulation and hemostasis; or one or more specific maikers of cardiac 
injury, one or more markers related to blood pressure regulation, and one or more madcers 
related to coagulation and honostasis, where each of these particularly prefared panels 
may optionally comprise one or more non-specific markers of tissue injury, markers 
related to inflammation, mariners related to pulmonary injury, and/or markers related to 
apoptosis. 

[0153] In similar fashion, a panel may comprise a plurality of maikera selected to 
diagnose, and/or distinguish amongst a plurality of, cerebrovascular disorders. In these 
aspects related to cerebrovascular disease, prefened marker panels comprise a phirality of 
maikers independently selected £pom the groiq> consisting of specific maikars of neural 
tissue injury, markers related to blood pressure regulation, maikers related to coagulation 
and hemostasis, matkers related to inflammation, and markers related to apoptosis. 
Preferably, such a panel comprises markerd ftom two, toee, four, or five different 
members of this gioi^. 

[0154] The following table provides an exemplary list of markers for use in the 
methods described herem: 



Marker 
Myoglobin 
• E-selectin 
VEGF 

Troponin I and complexes 
Troponin T and complexes 

Annexin V 

B-enolase 



Classification 
Nonspecific tissue iajury 
Nonspecific tissue injury 
Nonspecific tissue injury 
Myocardial injury 
Myocardial injury 
Myocardial injury 
Myocardial iajury 
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CK-MB 


Myocardial injury 


Glycogen phospharylase-BB 


M^cardialinjmy 


Heart type &tty add bindrng protein 


Myocardial injury 


Phosphoglycmc add mutase 


Myocardial injury 


S-lOOao 


Myocardial injury 


ANP 


Blood pressure regulation 


CNP 


Blood pressure regulation 


urotensin II 


Blood pressure regulation 


BMP 


Blood pressure regulation 


caldtonin gene related p^tide 


Blood pressure regulation 


arg-Vasopressin 


Blood pressure regulation 


Endothelin-1 


Blood pressure regulation 


Endalibelin-2 


Blood pressure regulaction 


Eiidotheliii-31 


Blood pressure regulation 


procaldtonjn 


Blood pressure regulation 


calcyphosiiiB 


Blood pressure regulation 


adreiumedullin 


Blood pressure regulation 


aldosterone 


Blood pressure regulation 


angiotensin 1 


Blood pressure regulation 


£Uigiotensin 2 


Blood pressure regulation 


angiotensin 3 


Blood pressure regulation 


Bradykinin 


Blood pressure regulation 


caldtonin 


Blood pressure regulation 


£ikbthdii>-2 


Blood pressure regulatiott 


Bidotheiljn-3 ' 


Blood pressure regulation 


Baiin 


Blood pressure regolation 


Urocnlatin 


Blood pressiure regulation 


Plasmin 




Thrombin 


Coagulation and hemostasis 


AntithroMbin-in 


Coagulation and hemostasis 


Fibrinogen 


Coagulation and hemostasis 


von Willebrand factor 


Coagulation and hemostasis 


D-dimer 


Coa^gulation and hemostasis 


PAI-1 


Coagulation and hemostasis 


PROTEIN C 


Coagulation and hemostasis 


TAFI 


Coagulation and hemostasis 


Fibrinopeptide A 


Coagulation and hemostasis 


Plasmin alpha 2 antiplasrain complex 


Coagulation and hemostasis 


Platelet factor 4 


Coagulation and hemostasis 


Platelet-derived growth factor 


Coagulation and hemostasis 


P-selectm 


Coagulation and hemostasis 


Prothrombin fragment 1+2 


Coagulation and hemostasis 


B-thrombogJobulin 


Coagulation and hemostasis 


Thrombin antithrombin III complex 


Coagulation and hemostasis 


Tbrombomoduhn 


Coagulation and hemostasis 


Thrombus Frecuisor Protein 


Coagulation and hemostasis 


Tissue &ctor 


Coagulation and hemostasis 
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basic c&lpoiiiii 1 


Vascular tissue 




Vascular tissue 




Vascular tissue 


CSKP2 


Vascular tissue 




Vascular tissue 


Fibniliul 


Vascular tissue 


LTBP4 


Vascular tissue 


stnooA muscle myosin 


VaflfMilar fissile 


transgelin 


Vascular tissue 


Carboxjftfinuiiial propeptide of type 




I procollsgcn (PICP) 


Collagen synthesis 


Coll&gco. cflii}oxyt63iiiiii3l telopcptidfi 




(ICTP) 


Collagen degradatioi 


Olutathioooe S TFansferase 


InflamTnatmry 


HIFl ALPHA 


Inilammatoiy 


ILrlO 


Ijoflammatory 


IL-l-Beta 


TnHammatntry 


IL-lra 


IhflaiDmatoiry 


IL-6 


Tnflammfltnry 


IL-8 


inflammatory 


Lysophospliatidic Ecid 


inflammatory 


MDA-modified LDL 


■ Liflanunatoty 


Human neutrophil elastase 


inflammatory 


C-reactiv6 protetn 


inflammatory 


Insulin-like growtii factor 


hiflammatory 


Inducible nitric oxide syn^ase 


Inflammatory 


intracellular adhesion molecule 


inflammatory 


Lactate dehydrogenase 


inflammatory 


MCP-1 


inflammatory 


MDA-LDL 


inflammatory 


MMP-1 


fiiflammatory 


MMP-2 




MMP-3 


Inflammatory 


MMP--9 


Inflammatory 


TIMP-1 


Mlammatory 


TIMP-2 


inflammatory 




Inflammatory 


n-acetj^ aspartate 


Inflammatory 


HTNTP "Qfv^ffmtnr QiTn*rfotniKr X^ArrkKiar 1 A 

Lrtr Rnwpvx oupcriiiuiuy jyicmDcr lA 


Jnflammatory 


Tbransfoi iuin§ growflh &ctor hota 




Tumor necrosis &ctor alpha 


Inflammato^f 


Vascular cell adhesion molecule 


inflammatory 


Vascular endothelial growth factor 


Inflammatory 


cystatin C 


inflammatory 


substance P 


Inflammatory 


Myeloperoxidase (MPO) 


Inflanomatory 


macrophage inhibitory &ctor 


Inflammatory 


Fibronectin 


Inflammatory 
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cardiotropbin 1 


Inflammatory 


Haptoglobin 


Ttiflammatnyy 


PAPPA 


laflammatory 


s-CD401igaiid* 


Inflammatory 


HMG 


Inflammatory 


IL-1 


inflammatory 


IL-2 


Inflammatory 


IL— 4 




IL-6 


Inflammatory 


IL-8 


inflammatory 


IL~10 


Inflammatory 


IL-U 


Tnflnmmatory 


IL-13 


Inflammatory 


IL-18 


Xnfiammstory 


Eosinophil catiomc protein 
Mast cell tryptase 


Inflammatory 
Inflammatory 


VCAM 


InfUunmatory 


sICAM-1 


Inflammatory 


TNFa 


Inflammatory 


Osteoprotegerin 


Inflammatory 


Prostaglaadin D-synthase 


Inflammatory 


Prostaglandin E2 


inflammatory 


RANKligand 


Inflammatory 


HSa'-60 


XnflamoEnatory 


Serum Amyloid A 


Inflammatory 


s-iL 18 receptor 


Inflammatory 


S*"iXj~I receptor 




S-TNFP55 


Inflammatory 


S-TNFP75 


Inflammatory 


TGF-beta 


Inflammatory 


MMP-11 


Inflammatory 


BetaNGF 


Inflammatory 


CD44 
EGF 


Inflammatory 
Inflammatory 


&-selectin 


Inflammatory 


Fibrondctin 


Inflammatory 


Neutrophil elastase 


Pulmonary injury 


KL-6 


Pulmonary injury 


LAMPS 


Pulmonary injury 


LAMPS 


Pulmonary injury 


Lung Surfactant protein A 


Pulmonary injury 


Lung Swfactant protein B 


Pulmonary injury 


Lung SnijEactait protein C 


Pulmonary ii^uiy 


Lung Surfactant protein D 


Pulmonary iiqury 


phospholipaseD 


Pulmonary injury 


PLA2G5 


Pulm^onaiy injury 


SFTPC 


Puhnonaiy injury 
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MAPKIO 
KCNK4 
KCNK9 
KCNQ5 
14-3-3 
4.1B 
AP0E4.1 
myelin basic piotdm 
Atcophin 1 
brain Derived nmtotrophic fitctor 
Brain Fatty add binding protem 
brain tubulin 
CACNAIA 
Calbindin D 
Calbrain 
Carbonic aohydrase XI 
CBLNl 
Cerebellinl 
Chimeiin 1 
Chiin6rin2 
CHNl 
CaEEN2 

Gliaiy nemobophic factor 
CK-BB 
CRHRl 
Otau 
DKPIA 
GFAP 
GFM6B 
GPR7 
GPRS 
GRIN2C 
GRM7 
HAPIP 
HIP2 
LDH 
Myelin basic protein 
NCAM 
NT-3 
NDPKA 
Neural cell adhesion naolecule 
NEUR0D2 
Neurofilimcnt L 
Neuroglobin 
neuromodulin 
Neuron specific enolase 
N6aiDpq)tide Y 



Neural tissue injury 
Neural tissue iqury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neund tissue injury 
Neural tissue injury 
Nepal tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue iigury 
Neural tissue injury 
Neontl tissue injury 
Neural tissue injury 
Neural ti^ue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Nairal tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
Neural tissue injury 
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Neurotensin 


Neural tis 


sue injury 


NeHTotrophin 1^^,4 


Neuraltis 


sue injury 


NRG2 


Neural tis. 


sue injury 


PACE4 


Neural tis 


sue injury 


phosphoglycerate mutase 


Neural tis 


sue injury 


PKC gamma 


Neural tia 


sue injury 


proteolipid protdn 


Neural tisj 


sue injury 


PTEN 


Neural tis 


sue injury 


PTPRZl 


Neural tis 


sue injury 


RGS9 


Neural tisj 


sue injury 


RNA Binding protein Regulatory Subumt 


Neural tisj 


sue injury 


S-100^ 


Neural tissue injury 


SCA7 


Neural tissue injury 


secretagogin 


Neural tissue injury 


SLC1A3 


Neural tis 


sue injury 


SORLl 


Neural tissue injury 


SREB3 


Neural tis 


>ue injury 


STAC 


Neural tis 


>ue injury 


STXIA 


Neural tissue injury 


STXBPl 


Neural tis£ 


me injury 


Syntaxin 


Neaural tissue injury 


thrombomodulin 


Neural tissue injury 


transthyretin 


Neural tissue injury 


adoiylate kinase- 1 


Neural tisaieinjnry 


BDNF* 


Neural tias 


reinjury 


neurokinin A 


Neural tissue injury 


s-acetyl Glutathione 


£^optosis 


cytochrome C 


apoptosis 


Caspase3 


apoptosis 


CathepsinD 


apoptosis 


o-spectrin 


^optosis 



[0155] Assay Measurement Strategies 

[01561 NurncTOUs methods and devices are weU known to the skilled artisan &T the 
detection and analysis of the markers of the instant invention. Witibi regard to polypeptides 
or proteins in patient test samples, immunoassay devices and methods are oftoi used. See, 
e.g., U.S. Patente 6,143,576; 6,113,855; 6,019,944; 5,985,579; 5,947,124; 5,939,272; 
5,922,615; 5,885,527; 5,851,776; 5,824,799; 5,679,526; 5,525,524; and 5,480,792, each of 
whidi is hereby incorporated by reference in its eaitirety, includmg all tables, figures and 
clainns. These devices and methods can utilize labeled molecules in various sandwich, 
competitive, or non-competitive assay formats, to generate a signal that is related to the 
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presence or amoimt of an analyte of interest. Additionally, certain methods and devices, 
such as biosensors and optical immunoassays, may be employed to determine ttie presence 
or amount of analytes without the need for a labeled molecule. See, e.g., U.S. Patente 
5,631,171; and 5,955,377, each of ^ch is heretoy incorporated by refeorence in its 
entirety, including all tables, figures and claims. One skilled in the art also recognizes HasX 
robotic instrumentation including but not limited to Beckman Access, Abbott AxSym, 
Roche ElecSys, Dade Bdning Stratus systems are among the immunoassay analyzers that 
are enable of perfonning the immunoassays taught herein. 

[0157] Preferably the markers are analyzed using an immunoassay, although o&er 
methods are well known to those skilled in the art (for example, the measurement of 
marker RNA levels). The presaice or amount of a marker is generally determined using 
antibodies specific &r each marker and detecting specific binding. Any suitable 
immunoassay may be utilized, for example, enzyme-linked immunoassays (ELDS A), 
radioimmunoassays (RIAs), competitive binding assays, and the like. Specific 
immunological binding of the antibody to the marker can be detected directly or indirectly. 
Direct labels include fluorescent or luminesceat tags, metals, dyes, radionuclides, and the 
like, attached to the antibody. Indirect labels include various enzymes well known in the 
art, sudi as alkaline phosphatase, horseradish peroxidase and the like. 

[0158] The use of immobilized antibodies specific for the markers is also 
contonplated by the present invention. The antibodies could be immobilized onto a 
variety of solid supports, such as magnetic or chromatographic matrix particles, the 
surface of an assay place (such as microtiter wells), pieces of a solid substrate material or 
membrane (such as plastic, nylon, paper), and the like. An assay strip could be prepared 
by coating the antibody or a plurality of antibodies in an array on solid support. This str^ 
could then be dipped into the test sample and ttien processed quickly through wa^es and 
detection stqjs to generate a measurable signal, such as a colored spot 

[0159] The analysis of a plurality of markers may be carried out separately or 
simultaneously with one test sample. For separate or sequential assay of markers, suitable 
apparatuses include clinical laboratory analyzes such as the ElecSys (Roche), fee AxSym 
(Abbott), tiie Access (Beckmaa), the ADVIA® CENTAUR® (Bayer) itnmunoassay 
systems, tiie NICHOLS ADVANTAGE® (Nichols Institute) immunoassay system, etc. 
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Preferred apparatuses or protein ciaps perform simultaaeoos assays of a plurality of 
markers on a smgle surfece. Particularly useful physical fimnats comprise surfeces having 
aplurality of discrete, addressable locations for the detection of aplurality of difierent 
aualytes. Such formats include protein miaxjanays, or **protein chips" (see, e.g., Ng and 
Hag, J. Cell Mol. Med. 6:329-340 (2002)) and certain capillary devices (see, e.g„ U.S. 
Patent No. 6,019,944). In these embodiments, each discrete surface location may 
comprise antibodies to immobilize one or moie analyte(s) (ag., a marker) for detection at 
each locatioa Surfecra may alternatively comprise one or more discrete particles (e.g.. 
miciopardcles or nanopaiticles) immobilized at discrete locations of a surface, where the 
microparticles comprise antibodies to irmnobilize one analyte (e.g., a marker) for 
detection. 

[01601 Several markers may be combined into one test for effident processing of a 
multiple of samples. £a addition, one stalled in the art would irecogoize the value of 
testing multiple samples (for example, at successive time points) fiom the same individual 
Such testing of serial samples will allow the identification of changes in maiker levels 
overtime. Increases or decreases m marker levels, as well as the aibseaice of change in 
maiker levels, would provide useful information about the disease status that includes, but 
is not hmited to identifying the approximate time from onset of the event, the presence and 
amount of salvagable tissue, the appropriateness of drug therapies, the effectiveness of 
various therapies as indicated by reperfiision or resolution of symptoms, differentiation of 
the various types of ACS, identification of the severity of the event, identification of the 
disease severity, and identification of the patient's outcome, including risk of future 
events. 

[G161] A panel consisting of the markers referenced above may be constructed to 
provide relevant infbrmatiwi related to differential diagnosis and/or prognosis. Such a 
pauel may be constucted using 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, or more or individual 
markers. The analysis of a single marker or subsets of markers comprising a larger panel 
of markers could be carried out by one skilled in the art to optimize clinical seaisitivity or 
specificity ia various clinical settings. These include, but are not limited to ambulatory, 
urgent care, critical care, intensive care, monitoring unit, inpatient, outpatient, physician 
oftice, medical clinic, and health screening settings. FurthermoFe, one skilled in the art 
can use a singje maiker or a subset of maikears conq>rising a larger panel of maikers in 
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combination with an adjustment of the diagtK>stic threshold in each of the aforementioned 
settings to optimize climcal sensitivity and specificity. The clinical sensitivity of an assay 
is defined as the percentage of those with the disease that the assay correctly predicts, and 
tiie specificity of an assay is defined as the perc^tage of fiiose without the disease that the 
assay correctly predicts (Tietz Textbook of Climcal Chemistry, 2*^ edition, Cail Burtis and 
Edward Ashwood eds., W.B. Saunders and Company, p. 496). 

[0162] The analysis of markers could be carried out in a variety of phj^ical formats as 
well. For exaaq>le, the use of microtiter plates or automation could be used to facilitate 
the processing of large numbers of test samples. Alternatively, single sample formats 
could be developed to &ciUtate immediate treatment and diagnosis in a timely i&sbion, for 
example, in ambulatory transport or emergency room settings. 

[0163] In another embodiment, the presoit invention provides a kit for the analysis of 
markers. Such a kit preferably comprises devises and reagents for the analysis of at least 
one test sample and instructions for performing the assay. Optionally the kits may contain 
one or more means for using information obtained fixjm mnrnmoassays performed for a 
marker panel to rule in or out certain diagnoses. 

[0164] Selection of Antibodies 

[0165] The generation and selection of antibodies may be acconqplished several ways. 
For example, one way is to pnrijEy polypeptides of interest ai to synthesize the 
polypeptides of interest using, e.g., solid phase peptide synfliesis methods well known in 
the art. See, e.g., Guide to Protein Piafication, Murray P. Deutcher, od, Metk EnzymoL 
Vol 182 (i 990); SoUd Phase Peptide Synthesis, CSieg B. Fields ed., Metk EnzymoL Vol 
289 (1997); Kiso et al. Chem. Pharm. Bull. (Tokyo) 38: 1 192-99, 1990; Mostafavi et al. 
Biomed. Pept. Proteins Nucleic Adds 1: 255-60, 1995; Fujiwara a/., Chan. Pharm. 
Bull. (Tokyo) 44: 1326-31, 1996. The selected polypeptides may then be injected, for 
example, into mice or rabbits, to generate polyclonal or monoclonal antibodies. One 
skilled in the art will recognize that many procedures are available for the production of 
antibodies, for example, as described in Antibodies, A Laboratory Manual, Ed Harlow and 
David Lane, Cold Spring Harbor Laboratory (1988), Cold Spring Harbor, N.Y. One 
skilled in the art will also appreciate that binding fiagments or Fab fi^agments which mimic 
antibodies can also be pr^ared fiom genetic information by various procedures (AntiTwdy 
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Bogineering: A Practical Approach (Borrebaeck, C, ed), 1995, Oxford Univearaity Press, 
Oxford; J. Immunol. 149, 3914-3920 (1992)). 

[0166] In addition, numerous pitolications have reported the use of phage display 
technology to produce and screen libraries of polypeptides for binding to a selected target. 
See, e.g, Cwirla et al., Proc, Natl. Acad. Sci, USA 87, 6378-82, 1990; Devlin et al.. 
Science 249, 404-6, 1990, Scott and Smith, Science 249, 386-88, 1990; and Ladner et al., 
U.S. Pat. No. 5,571,698. A basic concept of phage display methods is the establishment of 
a physical association between DNA encoding a polypeptide to be screened and the 
polypeptide. This physical association is provided by the phage particle, which displaj^ a 
polypeptide as part of a capsid enclosing the phage genmne which encodes the 
polypeptide. The establishment of a physical association between polypeptides and Iheir 
genetic matedal allows simultaneous mass soreening of very large numbers of phage 
bearing different polypeptides. Phage displaying a polypeptide with afGnity to a target 
bind to the target and these phage are enriched by affinity screening to the target. The 
identity of polypeptides displayed from these phage can be determined from their 
respective genomes. Using these methods a polypeptide identified as having a binding 
affinity for a desired target can then be synthesized in bulk by conventional means. See, 
e.g., U.S. Patent No. 6,057,098, v/bicb. is hereby incorporated m its entirety, including all 
tables, figuresj, and claims. 

[0167] The antibodies that are generated by diese methods may then be selected by 
first screening for affinity and specificity with Ihe purified polypeptide of int^st and, if 
required, comparing the results to the affinity and specificity of the antibodies with 
polypeptides that are desired to be excluded fi:om binding. The screening procedure can 
involve immobilization of the purified polypeptides in sq)arate wells of microtiter plates. 
The solution containing a potential antibody or groups of antibodies is then placed into the 
respective microtiter wells and incubated for about 30 min to 2 h. The microtiter wells are 
then washed and a labeled secondary antibody (for example, an anti-mouse antibody 
conjugated to alkaline phosphatase if the raised antibodies are mouse antibodies) is added 
to the wells and incubated for about 30 min and thcai washed. Substrate is added to the 
wells and a color reaction will appear whets antibody to the immobilized polypeptide(s) 
are pres^t^ 
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[0168] The antibodies so identified may thea be further analyzed for afHmty and 
specificity in the assay design selected. la the development of immunoassays for a target 
protein, the purified target protein acts as a standard with which to judge the sensitivity 
and specificity of the immunoassay using the antibodies tibat have been selected. Because 
the binding affinity of various antibodies may differ certsin antibody pairs {e.g., in 
sandwich assays) may inteifere with one another sterically, etc., assay performance of an 
antibody may be a more important measure than absohite afOnity and specificity of an 
antibody. 

[0169] Those skilled in the art will recognize that many approaches can be tak^ in 
producing antibodies or binding fiagments and screening and selecting for affinity and 
specificity for the various polypeptides, but these approaches do not change the scope of 
the invention. 

[0170] Selecting a Treatment Refflrpep 

[0171] The appropriate treatments for various types of vascular disease may be large 
and divarse. However, once a diagnosis is obtained, the cliniciaQ can readily select a 
treatment regimen that is con:q)atible with the diagnosis. Accordingly, the present 
invention provides methods of eariy differential diagnosis to allow for appropriate 
intervention in acute time windows. The skilled artisan is aware of ^ropriate treatments 
for numerous diseases discussed in relation to the methods of diagnosis described herein. 
See. e.g.. Merck Manual of Diagnosis and Uiere^, 17* Ed. Merck Research Laboratories, 
WMtehouse Station, NJ, 1999. 

[0172] The &>llowing provides a brief discussion of ^ditional exemplary markers for 
use in identifying suitable marker panels by the methods d^cribed herein. 

[0173] Examples: 

[0174] Example 1: Selection of Maikars for a Stroke Panel . A set of samples &om 
patients diagnosed with stroke and normal healthy donors were assayed for several 
markers of potential utility. No individual marker has sufficient clinical utility to diagnose 
stroke. The methods described above were used to determine the optimum markers for use 
in a panel of maricors. The data was separated into diseased and non-diseased groups. The 
indicator fhncttons were selected to be ramp ibnctions for all markers. The objective 
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function was chosen to be the product of the area, the koee, the specificity at 92.5% 
sensitivity and the sensitivity at 92.5% specificity. Hie initial simplex was randomly 
distributed about a vertesx derived fi:om the univariate analysis. Using tiw downhill simplex 
method with simulated annealing a local minimum was found that masdinized the 
objective function. For contribution for each analyte was calculated by setting the 
weighting parameter to zero and calculating the change in the objective function. This 
process was repeated 50 times. The mafkers were tanked by their average contribution 
over the 50 optimizations. The ROC curves for the initial vertex and an optimization are 
shown in figure 10. The ranking of the marker contributions is ishown in figure 1 1 . The 
lowest half of the markers were removed from the panel and the process was repeated. 
Figure 12 and 13 show the same mformation as in figures 10 and 1 1 but for the 19 marker 
panel. The lowest 9 markers were removed from the panel and the process was repeated. 
Figures 14 and 15 show the same mformation as in Figures 10 and 1 1 but for the 10 
marker panel. The lowest 5 markers were removed from the panel and the process was 
repeated a final time. Figure 16 and 17 show the same information as in figures 10 and 1 1 
but for the 5 marker panel. The individual ROC curves of the final 5 markers are shown m 
figure 9, The order of tiie contoibution does not match the order of the area of the 
individual ROC curves. A marker wifh poorer univariate utility may have greater utility 
when used in a panel. The area of the ROC curve deceases with decreasing panel size. 

[0175] Eyam ple 2: Tmprovement in diagnosis of AMI Utilizing Changes in Marker 
Levels . Data from a clinical study &om patients presenting with chest pain with serial 
draws from each patient was analyzed using the methods described in this document. The 
data was first analyzed without using daived markers. The data was again analyzed 
utilizing derived markers tiiat were related to tiie change in nuuk^ value firom the initial 
value. The ROC curves firom both optimized panel responses are shown in figure 18. The 
data cleariy illustrates the utility of the change in markeis to improve the diagnostic ability 
of panels in acute disease states. The method was also applied to detennine the best 3 and 
2 marker panels, and the results are also shown m figure 18. Figure 19 shows the 
contributions of tiie six AMI markers. Myoglobin, vi^le not a specific marker for AMI is 
a small molecule and the first marker of the three to elevate after AML Tnl is a specific 
marker for AMI, but is released more slowly. The method was not aware of this but still 
chose Tnl value arui change in Myoglobin. 
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[01761 Example 3: Simultaaeous Optimization of Two Criteria . Ie this example 
known stroke samples are analyzed with both stroke mimics and NHD samples in the non- 
diseased set. There are about 50 mimics and about 500 NHD samples, so the weighting is 
heavily in tsvox of optimizing results for KEID san^Ies. Aftsx optimization the panel 
response is fqjplied to a test set stroke vs. mimics and stroke vs. NHD. Similarly the data 
was optimized on stroke vs mimics, and the panel response was ^lied to as test set 
stroke vs. NHD and stroke vs NHD and mimics. Table 2 shows Hie average results of 
sanrple runs applied to the optimization sets and liien to the test sets. The effectiveness of 
the test is poor with respect to mimics. Two more optimizatioins were made as before, but 
fbis time a second group of data is sinmltaneously optimized. The second group consists of 
the stroke saaq>le8 and the mimics. Table 2 also shows the average results of sample runs 
^plied to the optimization set and vifbea the panel response is appUed to the two test sets. 
The efifectiveness of llie test witii respect to mimics is now improved. 

[OITT] While the invention has been described and exemplified in sufficient detail for 
tibiose skilled in this art to make and use it, various alternatives, modifications, and 
improvements should be apparent without departing fixnn the spirit and scope of the 
invention. 

[0178] One skilled in the art readily appreciates that the present invention is well 
adapted to carry out the objects and obtain tiie ends and advantages mentione4 as well as 
those inherent therein. The examples provided herein are representative of preferred 
embodiments, are exemplary, and are not intended as limitations on the scope of the 
invention. Modifications thjetdn and other uses will occur to those skilled in the art. These 
modifications are encompassed within Ihe spirit of the invention and are defined by Ihe 
scope of the claims. 

[0179] It will be readily apparent to a person skilled in the art that varying 
substitutions and modifications may be made to the invention disclosed herein without 
departing from the scope and spirit of the invention. 

[0180] All patents and publications mentioned in the specification are indicative of the 
levels of those of ordinary skill in the art to which the invention pertains. All patents and 
publications are herein incorporated by reference to the same extent as if each individual 
publication was specifically and individually indicated to be incorporated by refereoce. 
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10181] The invention illustratively descsribed herein suitably may be practiced in the 
absence of any elemeot or elements^ limitation or limitations which is not specifically 
disclosed herein. Thus, for exsanple, in each instance herem any of the teims 
"comprising", "consisting essentially of and "consisting of may be replaced with either 
of the other two terms. The temis and expressions which have beeai onployed ate used as 
terms of description and not of limitation, and there is no intention that in the use of such 
terms and e!q)ressions of excluding any equivalents of the features shown and described or 
portions thereof, but it is recognized that various modifications ate possible within the 
scope of ttie inv^tion claimed. "ITnis, it should be undesrstood that althougji the present 
invention has been specifically disclosed by preferred embodiments and optional features, 
modification and variation of the concq)ts herein disclosed may be resorted to by those 
skilled in the art, and that such modifications and variations are considered to be within 
Hhe scqpe of this invention as defined by Has attended claims. 

[0182] Othea: embodiments are set forth within the following claims. 
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We claim.: 



1. A method of identifyiiig a panel of markacs for diagnosis or prognosis of a disease 
or a condition, conqmsing: 

a) caloDlating a panel response for each patient in a set of diseased patients 
and in a set of non-diseased patients, said panel response being a function of value of each 
of a plurality of markers in a panel of maricers; 

b) calculating a value for an objective function, said objective fiinction being 
indicative of an effectiveness of the panel; and 

c) iterating steps a) and b) by varying at least one of parameters relating to 
said panel response fonctian and a sense of each madcer to facilitate optimization of said 
objective function. 

2. The method according to claim 1, whierein said objective functioo is a measure of 
aa overlap of panel res^nses of diseased patients and panel responses of non-diseased 
patients. 

3. The method accordmg to ckhn 1, whexein said panel response is a fimctioin of 
value of an indicator for each of a plurality of markers in a panel of markers and a 
weightiiig coefficient for each marker, said indicator being a mapping, for each of said 
plurality of markers;, of marker levels, said mapping being accarding to an indicator 
function; and 

wherein said iterating includes varying at least one of said weighting 
coefficieaats, parameters rdating to said indicator fonction, and a sense of each marker to 
facilitate optimization of said objective function. 

4. The method according to claim 3, wherein each indicator has a first value for 
marker levels below a cutoff region and a second value for marker values above a cutoff 
region, said cutoff region being defined by a location and a length. 

5 . The method according to claim 4, wherein said parametea:^ include said location of 
said cutoff region and said length of said cutoff region. 
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6. The method accoidmg to claim 4, wherein said Ifengfli of said cutoff region is zeax). 

7. the method accoidiag to claim 4, wberdu said length of said cutoff region is 
greater than zero. 

8. The method aoooaxling to claim 7, herein said indicatorB have values between said 
first value and said second value for marker levels wilimi said cutoffregion. 

9. The method accoidifig to daim 8, wherein said indicators have valu^ varying 
linearly &om said first value to said second value across said cutoffregion. 

10. The mediod according to claim 8, wheron said indicatorB have values varying non- 
linearly fix>m said first value to said second value across said cutoff region. 

1 1 . The method according to claim 1 0, wherein said non-linear variation is indicative 
of an error fimction of a distribution of marker values of diseased patients and an error 
function of a disfribution of marker values of non-diseased patients within said cutoff 

region. 

12. The method according to claim 3, wherein said calculating a panel response 
includes calculating, for each patient, Dwil^, where w is a weigjiting coefficient for a 
maimer i, I is the indicator value for the marker i, and X is a summation over all of said 
plurality of markers. 

1 3. The method accoFdis^ to claim 1 , wherein said calculating a value for an objective 
fimction includes generating a receiver operating characteristic (ROC) curve for said panel 
response, said ROC curve bemg indicative of a sensitivity of said panel respcaise as a 
function of one minus a specificity of said panel response. 

14. The method according to claim 13, wherein said objective function is associated 
vnth an area under said ROC curve. 
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1 5. The mdhod according to claim 13, wherdn said objective function Is associated 
wifli a knee of said ROC curve. 

16. The method accordiiig to claim 13, whetein said objective function is associated 
witti a seiisitivity at a selected specificity level 

17. The method according to claim 13, whoreia said objective fiinctioii is associated 
with a specificity at a selected sensitivity level. 

1 8. The method according to claim 13, wherein said objective ftinction is associated 
with two or more of an area under said ROC curve, a knee of said ROC curve, a sensitivity 
at a selected specificity level, and a specificity at a selected sensitivity level. 

19. The method according to claim 13, wherein said iterating constrains at least one of 
an area mder said ROC curve, a knee of said ROC curve^ a sensitivity at a selected 
specificity levd, and a specificity at a selected sensitivity level to be above about 0.9. 

20. The method according to claim 1 , further compiising: 

d) removing at least one of said maikers fi:Dm said panel; 

e) calculating a vahie of said objective fimction; and 

i) ddmnining a contribution of said at least one of said maikeocs to said 
objective function based on a result of step e). 

21. The method according to claim 20, jEurther comprising: 

g) repeating steps d) throng f) by removing a difBareait at least one of said 
markers from said panel; aod 

h) eliminating a marker from said panel of matkets in accordance with said 
contribution of said marker to said objective function. 

22. The method according to claim 1 , fiirther comprising; 

d) removing at least one of said markers from said panel; 

e) itearating steps a) and b) by varymg parameters relating to said panel 
response function to fecilitate optimization of said objective fimction; and 
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f) detennming a contribution of said at least one of said loatkers to said 

objective function based on a result of step e). 

23 . The method according to claim 22, further comprising: 

g) i^)eatiiigste[»d)1hroughQ by removing a dififexenft at le^ 
maikers bom said panel; and 

h) eliminating a marker fiom said panel of markers in aocoidance with said 
contribution of said marker to said objective fimctioa 

24. A system for identifying a panel of maikers for diagnosis or prognosis of a disease 
or a conditiQn, compdsing: 

means for calculating a panel response for each patient in a set of diseased patients 
and in a set of non-diseased patioats, said panel response being a fonction of value of each 
of a plurality of markers in a panel of markears; 

means for calculating a value for an objective function, said objective function 
being indicative of an effectiveness of said panel; and 

means for iteratively activating said means for calculating a panel response and 
said means for calculating a value for an objective fknction, by varying at least one of 
parameters relating to said panel response fiinction and a sense of each marker to facilitate 
optimization of said objective fimction. 

25. The system according to claim 24, wherein said objective fimction is a measure of 
an overlap of panel responses of diseased patients and panel responses of non-diseased 
patients. 

26. The system according to claim 24, wherein said panel response is a fimction of 
value of an indicator for each of a plurality of markeis in a pmel of matke3:s and a 
weighting coefGLdent for each marker, said indicator being a mapping, for each of said 
plurality of markers, of maiker levels, said mitpping being according to an indicator 
fimction; and 

wherein said means for iteratively activating is adapted to vary at least one 
of said weighting coefficients, parameters relating to said indicator function, and a sense 
of each maiker to fecilitate optimization of said objective function. 
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27. The system according to claim 26, wherein each indicator has a first value for 
marker levels below a cutoff region and a second value for madcer values above a catoff 
region, said cutoff region being defined by a location and a lengtii. 

28. . The method accordiiig to claim 27, wherein said parameters include said location 
of said cutoff region and said length of said cutoff region. 

29. The system according to claim 27, wherein said length of said cutoff region is zero. 

30. The s>Btem according to claim 27, wherein said lengfli of said cutoff region is 
greater than zero. 

3 1 . The system according to claim 30, wherein said indicators have values between 
said first value and said second value for marker levels within said cutoff region, 

32. The system according to claim 31, wherein said indicators have values varying 
linearly fiom said first value to said second value across said cutoff region. 

33. The system according to claim 32, wherein said indicators have values varying 
non-linearly fiom said first value to said second value across said cutoff region. 

34. The system according to claim 33, wberem said non-linear variation is indicative 
of an error function of a distribution of marker values of diseased patients and an error 
fimction of a distcibution of marker values of non-diseased patients within said cutoff 
region. 

35. The sj^tem according to claim. 26, wherein said means for calculating a panel 
response is adapted to calculate, for each patient, Ew,!,, where w is a weighting coefficient 
for a marker i, I is the indicator value for the marker I, and H is a summation over all of 
said plurality of madcers. 
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36. The system according to claim 24, wherem said means for calculatiug a value for 
an objective fonction is ada^ited to generate a receiver operating characteristic (ROC) 
curve for said panel response, said ROC curve being indicative of a sensitivity of said 
panel resfponse as a fimction of one minus a specificity of said panel response. 

37. The system according to claim 36, wherein said objective function is associated 
with an area under said ROC curve. 

38. The system accocding to claim 36, wherein said objective function is associated 
witii a knee of said ROC curve. 

39. The system according to claim 38, wham said objective flmction is associated 
with a sensitivity at a selected specificity level. 

40. The system acooidiog to claim 36, wherein said objective fimction is associated 
witii a specificity at a selected sensitivity level. 

41. The system according to claim 36, wherein said objective function is associated 
with two or more of an area undar said ROC curve, a knee of said ROC curve, a sensitivity 
at a selected specificity level, and a specificity at a selected sensitivity level. 

42. The system according to claim 36, whaein said means for iteratively activating is 
adapted to constrain at least one of an area under said ROC curve, a knee of said ROC 
curve, a sensitivity at a selected specificity level, and a specificity at a selected sensitivity 
level to be above aibout 0.9. 

43 . The system according to claim 24, further oompiismg: 

means for detennining a contribution of said at least one of said markers to said 
objective fimction, said means for detennitting bedng adapted to remove at least one of 
said markers Scorn said panel aod to activate said means for calculating a vahie £br an 
objective function. 
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44, Hie system accordmg to claim 43, finrtheir comprismg: 

means for eliminatiiig a marker from said panel of markers in accoidance with said 
contribution of said maika: to said objective functiou, said means for eliminating beiag 
adapted to activate said means fat detennining a contributioii by removing a different at 
least one of said mariners &om said panel 

45, The system according to claim 24, ftnthea: comprising: 

means for determining a contribution of said at least one of said markers to said 
objective function, said means for detennining being adapted to remove at least one of 
said markers from said panel and to iteratively activate said means for calculating a panel 
response and said means for calculating a value for an objective function, by vaiying 
parameters relating to said panel response function to facilitate optimizatiQn of said 
objective function. 

46. The system according to claim 45, further coir^sing: 

means for eliminating a marker from said pand of markers in acooidance with said 
contribution of said marker to said objective function, said means for eliminating being 
adapted to activate said me^ for detenniidng a cantribntiom by removing a different at 
least one of said markers teem said panel 

47. A program product, comprismg rnachine readable program code for causing a 
machine to perform fbllowmg method steps: 

a) calculating a panel response Sar each patient in a set of diseased patiente 
and in a set of non-diseased patients, said panel response being a fhnction of value of each 
of a plurality of maikers in a panel of markeirs; 

b) calculating a value for an objective fimctian, said objective function being 
indicative of an effectiveness of said panel; and 

c) iterating steps a) and b) by varying at least one of parameters relating to 
said panel response function and a sense of each marker to f^ilitate optimization of said 
objective fimction. 
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48. The program product according to claim 47, wherein said objective fimctlon is a 
measure of an overlap of panel responses of diseased patients and panel responses of non- 
diseased patieats. 

49, The program product acco!rding to claim 47, wherein said panel response is a 
fimction of value of an indicator for each of a plurality of markers in a panel of markers 
and a weig^iting coefBcieait for each marker, said indicator being a mapping, far each of 
said plurality of marketrs, of marker levels, said mapping being according to an indicator 
function; and 

vdierein said iterating includes varying at least one of said weighting 
coefficients, parameters relating to said indicator fimction, and a sense of each marker to 
&cilitate optimization of said objective fimction. 

5 0. The program product accordiag to claim 49, wherein each indicator has a first 
value for marker levels below a cutoff regicai and a second value for marker values above 
a cutoff region, said cutoff region bdng defined by a location and a lengtli. 

5 1 . The program product according to claim 50, wherein said parameters include said 
location of said cutoff region and said length of said cutoff region. 

52. The program product according to claim 50, wherein said length of said cutoff 
region is zero. 

53 . The program product according to claim 50, wherein said length of said cutoff 
region is greater than zero. 

54. The program product according to claim 53, wherein said indicatois have vahies 
between said first value and said second value for raaikea: levels within said cutoff regtoo. 

55. The program product according to claim 54, wherein said indicators have values 
varying linearly Sxna. said first value to said second value across said cutoff region. 
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56. The program product accoidxog to claini 54, wherem said indicatois have values ' 
varying nan-linearly £rom said first value to said second value across said cutoff region. 

57. The program product according to claim 56, wherein said non-linear variation is 
indicative of m error fimction of a distribution of marker values of diseased patients and 
an error function of a distribution of marker values of non-diseased patients within said 
cutoff region. 

58. The pragram product according to claim 49, wberein said calculating a panel 
response includes calculating, for each patient, Ewjli, where w is a weighting coefficient 
for a marker i, I is the indicator value for the marker i, and Z is a summation over all of 
said plurality of markers. 

59. The program product aoooiding to claim 47, wherein said calculating a value for an 
objective fimction inchides generating a receiver opacating characteristic (ROC) curve for 
said panel response, said ROC curve being indicative of a sensitivity of said panel 
response as a fimction of one minus a specificity of said panel response. 

60. The program product accordiiig to claim 59, wherein said objective fmictioni^ 
associated with an area under said ROC curve. 

61 . The program product according to claim 59, wherdn said objective functiooa is 
associated with a knee of said ROC curve. 

62. The program product according to claim 59, wherein said objective function is 
associated with a sensitivity at a selected specificity level. 

63. The program product according to claim 59, whearein said objective function is 
associated witii a specificity at a selected sensitivity level. 

64. The program product according to claim 59, whearedn said objective function is 
associated with two or more of an area under said ROC curve, a knee of said ROC curve, a 
sensitivity at a selected specificity level, and a specificity at a selected sesisitivity level. 
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65. The program product according to claim 59, wherein said iterating constrains at 
least one of an area under said ROC curve, a knee of said ROC crave, a sensitivity at a 
selected spedfuxity level, and a specificity at a selected seositiyity level to be above about 
0.9. 

66. The program product according to claim 47, further comprising madune readable 
program code for causing a machine to per&mn. following method steps: 

d) rranoving at least one of said madceirs from said panel; 

e) calculating a value of said objective function; and 

f) determining a contribution of said at least one of said markers to said 
obj ective fimction based on. a result of step e). 

67. The program product according to claim 66, further comprising machine readable 
program code for causing a machine to perform following method steps: 

g) repeating steps d) through f) by removing a different at least one of said 
markers Scorn said panel; and 

h) eliminating a ma±^ fiom said panel of markers in accordance with said 
CQUttibution of said marker to said objective fonctioii. 

68. The program product according to claim 47, £uther contusing machine readable 
program code for caiuitig a machine to perform following meffaod steps: 

d) removing at least one of said markers from said panel; 

e) iterating steps a) and b) by varying parametears relating to said paiiel 
response fimction to facilitate optimization of said objective function; and 

f) determining a contribution of said at least one of said maikears to said 
objective function based on a result of step e). 

69. The program product according to claim 68, fiirthcrcoD]^ris^ 
program code for causing a machine to perform following method steps: 

g) repeating steps d) through f) by removing a different at least one of said 
markers from said panel; and 

h) eliminating a niaikerfix>m said panel ofmarkers in accordance with said 
contribution of said marker to said objective function. 
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70. The program product according to claim 47, wherein said machine readable code is 
embedded in a portable meto-. 

71 . The program product according to claim 70, wherein said portable meter is a 
fluorometra:. 

72. The program product according to claim 70, wherem said portable meter is a 
reflectometer. 

73 . The program product according to claim 47, wherein said machine readable code is 
embedded in a computer. 

1A. The progOTi product according to claim 73, wherem said computer is a portable 
computer. 

75. The program product according to claim 73, wherein said oon^uter is ada^ 
be accessed throu£^ a netwoik. 

76. The program product according to claim 75, wherein said network is the Ihtemet 

77. The prograin product according to claim 73, wherein said computer is adi^ted to 
be coiq)led to an analyzer. 

78. The program product according to claim 77, wherein said analyzer is an 
immunoassay analyzer. 

79. The program product according to claim 77, wherein said analyzer is a single 
nucleotide polymorphism detector. 

80. The program product according to claim 77, wherein said analyzer is adapted to 
sort and count similar and difTerent particles and cells. 
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81. A mdliod of identifying a panel of maikers for diagnosis or prognosis of a disease 
or a condition, comprising: 

a) selecting a panel of markers, said panel including a plurality of markers 
measured in a set of diseased patients and a set of non-diseased patients; 

b) defining a cutofifregiaa of marker levels fbr each of said plurality of 
matkets, said cutoff region having a location and alengfh; 

c) selectiiigawdghtingcoefEicieatfoTeaGhof saidplural^^ 

d) mt^ing, for each of said ptaralityofrnarke3S,niaiker levels torn 
indicator, each of said indicators lurving a first value for marker levels below said cutoff 
region and a second value for marker levels above said cutoff region; 

e) calculating a panel response for each patient in said sd of diseased patioots 
and in said set of non-diseased patients, said panel response being a fimction of value of 
said indicator for each marker and said weighting ooefQcient fbr each marker; 

f) calculating a value for an objective function, said obrjective function being 
indicative of an effectiveness of said panel; and 

g) iterating steps e) and f) by varying at least one of said location of said 
cutoff region, said length of said cutoff region, said wdghting coejeBcients, and a sense of 
each marker to facilitate optimization of said objective function. 

82. The method according to claim 81, wherein said objective function is a measure of 
an ovedap of panel responses of diseased patients and panel responses of non-diseased 
patients. 

83. The method according to claim 81, wherein said length of said cutoff region is 

84. The method according to claim 81, wherem said length of said cutoff region is 
glister than zero. 



85. The mettkod according to claim 84^ wherein said indicators have vahies between 
said first value and said second value for marker levels within said cuto£f region. 
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86. The method according to claim 85, wherein said indicators have values varying 
linearly from said first value to said second value across said cutoff fegion. 

87. The mediod according to claim 85, ^ordn said indicators have values varying 
non-linearly j&om said first value to said second value across said cutoff legion. 

88. The meftrad accoidmg to claim 87, wherein said non-linear variation, is indicative 
of an error function of a distribution of matter values of diseased patients and an error 
function of a distribution of marker values of non-diseased patients within said cutoff 
region. 

89. • The method according to claim 81, wherein said calculating a panel response 
includes calculating, for each patient, Swilj, where w is a weigjiting coefficient for a 
nmker i, lis the indicator value for the marker I, and £ is a summation over all of said 
phirality of markers. 

90. The method according to claim 8 1 , wfaeiein said calculating a value for an 
objective fimction includes generating a receiver operating chaxacteiistic (ROQ curve for 
said panel response; said ROC curve being indicative of a sensitivity of said panel 
response as a fonctioa of one minus a specificity of said panel response. 

91 . The method according to claim 90, wherein said objective fimction is associated 
with an area under said ROC curve. 

92. The method according to claim 90, wherein said objective fimction is associated 
with a knee of said ROC curve. 

93 . The method accotdmg to claim 90, wherem said objective fimction is associated 
with a sensitivity at a selected specificity level. 

94. The method accordiiig to claim 90, wherein said objective fimction is associated 
with a specificity at a selected sensitivity level. 
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95. The mctfhod accoiding to claim 90, Tudierein said objective function is associated 
witiitwo ormore of aa area nnder said ROC curve, aknee of said ROC curve, a sesnsitivity 
at a selected specificity level, and a specificity at a selected saisitivity level. 

96. The method accoidii]g to claim 81, further comprism 

h) setting said weighting coefficient of at least one of said maikets to 
proximately zero; 

i) calculatmg a vahieofsaid objective fimction with lemaiDi^ 
coefQcieats; and 

j) detemuning a contribution of said at least one of said markers to said 
objective function. 

97. The method according to claim 96, ftirther coinprising: 

k) repeating steps h) through j) by setting said weighting coeffident of a 
different at least one of said markers to approximately zsto; and 

1) eliminating a maikerfix>m said panel ofmarkers in accordance with said 
contribution of said marker to said objective function. 

98. A method of diagnosing a subject for a disease or condition, comprising: 

a) measuring a value of each of a plurality of maikars in said subj ect; 

b) calculating a panel response for said subject, said panel response being a 
function of values of a plurality of maikears; and 

c) diagnosing said disease oc condition m said patient when said panel 
response is greater than a predetermined cutoff . 

99. A program product, comprising machine readable program code for causing a 
machine to perfixon ibllowing method steps: 

a) calculating a value of each of a plurality of markers in a subject; 

b) calculating a panel r^ponse for said subject, said panel response being a 
function of values of a plurality of matkers; tmd 

c) diagnosing a disease or condition in said patient when said panel response 
is greater than a predetennined cutoff. 
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100. The program product according to claim 99, wherein said machine readable 
program code is onbedded in a portable meter. 

101. The program product according to claim 100, wherein said portable meter is a 
fluorometer. 

102. The program product according to claim 100, wherein said portable meter is a 
reflectomcfter. 

103. The program product according to claim 99, wherein said maflV i inft readable 
program code is embedded in a computer. 

104. The program product according to claim 103, wherein said carnputer is a portable 
computer. 

105 . The program product according to claim 103, wherein said computer is adapted to 
be accessed Ihrough a network. 

106. The program product according to claim 105, wherein said network is the Jhtemet. 

107. The program product according to claim 103, wherein said computer is adapted to 
be coiqpled to an analyser. 

108. The program product according to claim 1 07, wherein said analyzer is an 
immunoassay aaialyzer. 

1 09. The program product according to claim 1 07, wherein said analyzer is a single 
nucleotide polymorphism detector. 

110. The program product according to claim 107, wheidn said analyzer is adapted to 
sort and count smiilar and different particles and cells. 
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111. A method of identifying a panel of markers for diagnosis or prognosis of a disease 
or a condition, comprising: 

a) identifying a cutoff region for each of a plurality of markers, said cutoff 
region being substantially centered about an overlap region of marker vahies for a set of 
diseased patients and a set of non-diseased patients, said cutoff region having a location 
andalengtt^ 

b) detemdning an effectiveness value of each of said plurality of maikers in 
distinguishing said set of duieased patients fiom said set of non-diseased patients; and 

c) defining a panel response as a jKmctionofsaide£[ectiveQess value of each 
msAxx and a measured level of each marker. 

1 12. The method acc»>rding to claim 1 1 1, wherein said cutoff region has a length of 
zero. 

113. The method according to claim 111, wherein said cutoff region has a non-zero 
length. 

1 14. The method according to claim 1 1 1, wherein said efifoctivenfiss value of each 
maiker is represented by an area under a ROC curve. 

115. The method according to daim 3, wherdn said indicator function is monotonic 
with marker value. 

1 16. Hie method aocording to claim 115, wherein said indicator function is one of the 
gfxtxep confdsting of: a ranq) function, a step fimction, and a sigmoid function. 

117. The method according to claim 3, wherein said indicator function is ad^ted to 
localize a marker vahie. 

1 1 8. The method according to claim 1 17, vrtietein said indicator function is one of the 
group consisting of: a triangle, a square, and Gaussian. 
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119. The mediod according to claim 1 , ^iiiecdn at least one of said plurality of markers 
is a derived marker. 

120. The mdhod according to claim 119, wherein said derived marker is Ihe ratio of 
two other maricers. 

121. The method according to claim 1, wherran said iterating includes using a downhill 
simplex method. 

122. The mefliod according to claim 121, wherein said iterating further includes 
simulated annealing. 

123 . The mefliod according to claim 122, wherein said simulated annealing includes 
perfonning a statistically sufiScient number of optimizations to evaluate a most common 

solution. 

124. The melfaod according to claim I, wherein said optimizatian is adapted to provide 
a stable solution. 

125. The method according to claim 124, 'Vidierein said ad^tation inctudes varying the 
madcer values by a random percentage. 

126. The method according to claim 124, wherein said adaptation includes varying one 
or more parameters of said panel function. 

127. The mdhod according to claim 124, wherein said adaptation includes gotierating a 
seed simplex about a minimum. 

128. The method according to claim 124, wherein said adaptation includes increasing an 
annealing temperature until an achieved solution is not recovered. 
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129. The method according to claim 98, wherein said calculating a panel response 
includes usmg a panel response function and parameters relating to said panel response 
fimctioD, said panel response function and parameters being determined by: 

a) calculating a panel response for each patient in a set of diseased patients 
and in a set of non-diseased patients, said panel response bdng a function of value of each 
of a plurality of madcers in apanel of markers; 

b) calculating a value for an objective function, said objective function being 
indicative of an effectiveness of the panel; and 

c) itoiating 5tq>s a) and b) by varying at least one of parameters relating to 
said panel response functioii and a sense of each maifker to facilitate optinuzatioa of said 
objective fiinctioiL 

1 30. The method according to claim 129, wherem said panel of markers is adapted to 
diagnose two or more diseases or conditions. 

131. The method according to claim 130, wherein different panel response fjomctions 
and parameters are used for each of said diseases or conditions. 

132. The method accprding to claim 98, whecein said at least one of said plurality of 
markers is a derived marker. 

133. The method according to claim 132, whereui said derived marker is a ratio of two 
other markers. 

1 34. The method according to claim 132, wherein said derived marker is indicative of a 
change over time in a marker. 

IBS. The method according to claim 98, wheirein said calculatmg a panel response 
includes calculating, for each patient, Ewflj, where w is a weighting coefGcient for a 
marker i, lis an indicator value for the marker 1, and E is a summation over all of said 
plurality of maricQTS, said indicator value being a mapping, for each of said plurality of 
markers, of marker levels, said mapping being according to an indicator fimction. 
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136. Hie method according to claim 135, wherein said indicator function is moootonic 
with marker value. 

137. The method according to claim 136, wherein said mdicator function is one of the 
group consistiiig o£ a wnp function, a step function, and a sigmoid function. 

138. The method according to claim 135, wherein sadd indicator function is adapted to 
localize a marker value. 

139. The method according to claim 138, wherein said indicator fimction is one of tiie 
group consisting of: a trimgle, a square, and Gaussian. 

140. The metibK)d according to claim 20, wherdn said removing includes setting a 
weighting coefQcient of said at least one of said markers to approximately zero. 

141 . The system according to claim 26, wherein said indicator fimction is monotonic 
with marker value. 

142. The system according to claim 141 , wherein said mdicator function is one of the 
group consisting o£ a tanap function, a step function, and a sigmoid function. 

143. The system according to claim 26, wherein said indicator function is adapted to 
localize a marker value. 

144. The system according to claim 143, wherein said indicator i^mction is one of the 
group consisting of: a triangle, a square, and Gaussian. 

145 . The system according to claim 43, wherein said removing includes setting a 
weighting coefficient of said at 1^ one of said madcers to approximately zero. 

1 46. The program product according to claim 49, wherem said indicator function is 
monotonic with marker value. 
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1 47. The program product according to claim 49, wherein said indicator function is one 
of the group consisting of: a ramp fimctiaD, a step fimction, and a sigmoid function. 

148. The program product according to claim 49, wherdn said indicator function is 
adapted to localize a marlcer value. 

149. The progjmi product according to claim 49, wheiroin said indicator function is one 
of the group consisting of: a triangle, a square^ and Gaussian. 

150. The progrsun product according to claim 99, wherein said calculating a panel 
response includes using a panel response function and parameters relating to said panel 
response function, said panel response function and parameters being detennined by: 

a) calculating a panel response for each patient in a set of diseased patients 
and in a set of non-diseased patients, said panel response being a fimctioii of value of each 
of a plurality of markers in a panel of markers; 

b) calculating a value for an objective function, said objective function being 
indicative of an effectiveness of the panel; and 

c) iterating steps a) and b) by varying at least one of parameters relating to 
said panel response fimction and a sense of each marker to &cilitate optimization of said 
objective function. 

151. The program product according to claim 150, wherein said panel of maikeis is 
adE^ted to diagnose two or more diseases or conditions. 

152. The program product according to claim 151, wherein different panel response 
functions and parameters are used for each of said diseases or conditions. 

153. The program product according to claim 99, wherem said at least one of said 
pluraUty of markers is a derived marker. 

154. The program product according to claim 153, wherein said derived marker is a 
ratio of two other mackers. 
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155. The program product according to claim 153, herein said dedved marker is 
iitdicative of a change over time in a maiker. 

156. The program product according to claim 99, wherein said calculating a panel 
response includes calculating, for each patient, Swili, where w is a weighting coefBcient 
for a marker i, I is an indicator value for the marker i, and E is a summation over all of 
said plurality of madcers, said indicator value being a mapping for each of said plurality 
of markers, of marker levels, said mapping beitig according to an indicator fimction. 

157. The program product according to claim 156, wherein said indicator fimction is 
monotonic with marker value. 

158. The program product accordmg to claim 157, wherein said indicator fimction is 
one of the group consisting of: a ramp function, a step fimction, and a sigmoid fimction. 

159. The program product according to claim 156, herein said indicator fimction is 
adapted to localize a marker value. 

160. The program product according to claim 159, wherem said indicator fimction is 
one of the group consisting of: a triangle, a square, and Gaussian, 

161. Tbt method according to claim 1 19, wherein said derived marker is iiuHcative of 
the diange in another marker over time, 

1 62. The method according to claim 119, wherein said derived marker is indicative of 

the change in said panel response over time. 

163. The method according to claim 132, wherein said dmved marker is indicative of a 
change over time in said panel respanse. 

1 64. The program product according to claim 66, wherein said removing includes 
setting a weighting coefficient of said at least one of said markers to approximately zero. 
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165. The method according to claim 1 , wherem said optimization is adapted to 
simultaneously at least one of optimize and constrain a phuality of objective fijnctions 
calculated fipom a plurality of groins of data. 

1 66. The system accoiding to claim 24, whemn said means for iterativdy activating is 
ad^ted to simultaneously at least one of optimize and constrain a plurality of objective 
fmctioQS calculated fiom a plurality of groups of data. 

167. The program product according to claim 47, viiierein said optimizatiou is adapted 
to simultaneously at least one of optimize and constrain a plurality of objective functions 
calculated fiom a plurality of groups of data. 

168. The method according to claim 134, \diea:ein said calculating a panel response is 
adapted to calculate a panel response without said derived marker when said derived 
marker is not available and to calculate a panel response using said derived marker when 
said derived maikea: is available. 

1 69. The method according to claim 1 68, wherein said panel response without said 
derived marker and said panel response using said derived marker are calculated according 
to different panel response functions and parameters rel^oting to said panel response 
Amotions. 

170. The mdhod according to claim 163, wherein said calculating a panel response is 
ad^ted to calcukte a panel response wiliiout said derived marker when said derived 
marker is not available and to calculate a panel response using said derived markea: when 
said derived marker is available. 

171. The method according to claim 170, y^erdn said panel response without said 
derived marker and said panel response using said derived marker are calculated according 
to different pand response functions and parameters relating to said panel respcmse 
functions. 
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